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16500 INDUCED CLIMATE CHANGE 


KEY WORDS: Agricultural economics; Agricultural engineering; Carbon 
dioxide; Climatic changes; Environmental effects; Environmental research; 
Temperature effects; Temperature rise; Water requirements; Water shortage 


ABSTRACT: The unrestrained combustion fossil fuels possibly the most 
important environmental issue facing mankind. The concentration CO, the 
atmosphere may double within the next years, producing changes amounts 
rainfall, increasing temperatures and changing the agricultural potential. Temperature 
changes will greatest (about 8°C) degrees latitude. doubling CO, 
concentrations will increase temperatures near the equator less than 2°C latitude 
40° about 3.6°C. the latitudes 37° 47° rainfall projected decrease while 
water requirements are estimated increase possibly much percent. The 
potential for future agricultural production the United States will depend largely 
the availability additional water. Quadrupling CO, concentrations projected 
increase temperature more than 6°C 40° latitude. 


REFERENCE: Hargreaves, George (Research Engr., Dept. Agric. and Irrig. 
Engrg., Utah State Univ., Logan, Utah), “Water Requirements and Man-Induced 
Climate Change,” Journal the Irrigation and Drainage Division, ASCE, Vol. 107, 
No. IR3, Proc. Paper 16500, September, 1981, pp. 247-255 


16502 FURROW ADVANCE RATES UNDER SURGE FLOW SYSTEMS 


KEY WORDS: Automation; Flow rates; Furrow irrigation; Furrow systems; 
Hydraulic properties; Intakes; Irrigation design; Time factors; Water flow 


ABSTRACT: Two field tests were conducted study the effect cycling furrow 
inflows advance rates. Earlier investigations had demonstrated the concept and 
termed “surge flow” irrigation based the hydraulic behavior the cycled, 
pulsed, inflow along the furrow. Surge flow treatments were compared with continuous 
flow tests having the same time-averaged flow rate. All the tests indicated that surge 
flow produced significantly faster advance rates. The primary affects were noted during 
the water application newly formed furrows with lesser differences subsequent 
irrigations. Temporal and spatial variations furrow intake rates were substantially 
reduced the surge flow practices indicating major advantage the technique. 


REFERENCE: Bishop, Alvin (Prof. Emeritus, Dept. Agr. and Irrigation Engrg., 
Utah State Univ., Logan, Utah 84322), Walker, Wynn R., Allen, Niel L., and Poole, 
Gregory J., “Furrow Advance Rates under Surge Flow Systems,” Journal the 
Irrigation and Drainage Division, ASCE, Vol. 107, No. IR3, Proc. Paper 16502, 
September, 1981, pp. 257-264 
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KEY WORDS: Basins; Dimensional analysis; Distribution; Flow rate; 
Hydraulics; Irrigation engineering; Mathematical models; Open channel 
flow; Soil moisture; Surface irrigation 


ABSTRACT: The procedures are examined for obtaining reasonable estimates 
distribution uniformities for wide variety flow rates, length run, infiltration 
characteristics, and flow resistance for the design and management level basins. The 
dimensionless solution advance and recession level basins was extended show 
the distribution uniformities for wide variety conditions. This 
transformed into two representations distribution uniformity that are more useful for 
designing and managing level basins. One solution displays the effects soil moisture 
deficit and the necessary infiltration opportunity time distribution uniformity. The 
other displays the effects field length and flow rate distribution uniformity. 


REFERENCE: Clemmens, Albert (Research Hydr. Engr., United States Water 
Conservation Lab., 4331 East Broadway Road, Phoenix, Ariz. 85040), Strelkoff, 
Theodor, and Dedrick, Allen R., “Development Solutions for Level-Basin Design,” 
Journal the Irrigation and Drainage Division, ASCE, Vol. 107, No. IR3, Proc. 
Paper 16503, September, 1981, pp. 265-279 


16513 SIMPLIFIED METHOD FOR RAINFALL INTENSITIES 


KEY WORDS: Duration curves; Observation; Precipitation (meteorology); 
Precipitation rate; Rain and rainfall; Rainfall intensity; Regression analysis; 
Rods; Time factors 


ABSTRACT: Procedures are presented for estimating maximum rainfall amounts for 
hours with return period years from monthly and observation day 
measurements. Procedures are also presented for converting these values other 
rainfall depth- duration amounts. The depth-duration amount varies with the fourth 
root the time duration and the sixth root the number years the return 
period. Observation-day rainfall converted 24-hour values multiplying 1.13. 
Mean monthly rainfall and the maximum monthly values from 30-year series are 
used estimate 10-year return period, 24-hour duration amounts. The 10-year, 24- 
hour amounts are then used estimate other depth-duration values. The relationships 
are good except protected areas such for the lee mountains areas 
surrounded mountains. 


REFERENCE: Hargreaves, George (Research Engineer, Dept. Agri. and 
Irrigation Engrg. Utah State Univ., Logan, Utah.), “Simplified Method for Rainfall 
Intensities,” Journal the Irrigation and Drainage Division, ASCE, Vol. 107, No. 
IR3, Proc. Paper 16513, September, 1981. pp. 281-288 


16518 BORDER IRRIGATION SYSTEMS 


KEY WORDS: Crop production; Crop yield; Design; Economic analysis; 
Irrigation; Optimization; Simulation models; Systems engineering; Water 
requirements; Water resources 


ABSTRACT: Using surface irrigation hydraulics simulation model, relationships were 
developed between water requirement efficiency and the system design variables. 
crop production function was utilized relate crop yield the water requirement 
efficiency. Gross returns from the crop, and the costs water, labor, ditch 
construction and crop production were considered the optimization problem. Several 
system constraints were incorporated into the design process. The generalized 
geometric programming technique was applied the optimal design border and 
basin irrigation systems. The design variables were: the length run, inflow rate into 
the border, time inflow, number lengths run, width the border, and the 
number border widths the field. 


REFERENCE: Reddy, Mohan (Post-Doctoral Fellow, Dept. Civ. Engrg., 
Colorado State Univ., Fort Collins, Colo.), and Clyma, Wayne, “Optimal Design 
Border Irrigation Systems,” Journal the Irrigation and Drainage Division, ASCE, 
Vol. 107, No. IR3, Proc. Paper 16518, September, 1981, pp. 289-306 


16507 HYDRAULIC TESTING PLASTIC FILTER FABRICS 


KEY WORDS: Bacteria; Drainage; Fabrics; Filters; Ground water; 
Hydraulics; Hydraulic tests; Laboratory tests; Permeability; Pipes; Plastics; 
Sediments; Soils; Ultrahigh frequencies 


ABSTRACT: The hydraulic characteristics nonwoven filter filters were investigated 
using laboratory permeameter simulate soil/fabric subsurface drainage system. 
utilizing information discharge, pressure, permeability, and filter properties, 
conclusions were drawn concerning the correlation between the soil and the plastic 
filters. sands with large concentrations fine particulate, the migration the fines 
was minimal, and the soil controlled the hydraulic response the system regardless 
the filter fabric used. fine soil where the water velocity greatly reduced, bacterial 
activity within the soil can cause erroneous results. Some design criteria are developed 
for the use plastic filter fabrics. 


REFERENCE: Chen, Yung Hai (Assoc. Prof. Civ. Engrg., Colorado State Univ., 
Fort Collins, Colo. 80523), Simons, Daryl B., and Demery, Phillip M., 
Testing Plastic Filter Fabrics,” Journal the Irrigation and Drainage Division, 
ASCE, Vol. 107, No. Proc. Paper 16507, September, 1981, pp. 307-324 
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WATER REQUIREMENTS AND MAN-INDUCED 
CLIMATE CHANGE 


INTRODUCTION 


irrigation project planning, provision frequently made for dry and wet 
cycles well other extremes. Irrigation facilities, including storage dams, 
canals, and related structures are often designed for projected useful life 
years. Economic benefits are estimated largely from review present 
yields, prices, and water requirements. Planning and design can significantly 
improved reliable estimates can made probable future changes climatic 
well economic conditions. Until now, little attention has been given 
the magnitude immediacy possible man-induced climatic changes the 
rates which such changes can expected occur. 

Although the magnitude future climate change has not been accurately 
determined, enough now known indicate need for increasing research 
efforts and for evaluating new methodologies for monitoring the effects present 
and future climatic conditions food and other agricultural production. 

This paper based the assumption that increasing use fossil fuels 
will produce fairly rapid global temperature increase due the increase 
CO, the atmosphere. The theory that cooling cycle probable rejected. 

From current models effects that are projected due increased atmospheric 
carbon dioxide (models currently being used the U.S. Dept. Energy 
referenced later this paper) and assuming that the magnitude predicted 
temperature increases are approximately correct, attempt made this 
paper estimate ‘resulting changes that may occur water requirements and 
water management irrigated and dryland agriculture. Potential evapotranspira- 
tion, ETP, used index water requirements. Justification presented 
for using estimating equation based only temperature and radiation. 
Temperature increases predicted climate models are used together with ETP 
estimate approximate increase water requirements. 

Some proposed research topics relating water management are presented. 
The overall purpose the paper promote study, research, discussion, 

Engr., Dept. Agr. and Irrigation Engrg., Utah State Univ., Logan, Utah 
84322. 

Note.—Discussion open until February 1982. extend the closing date one month, 
written request must filed with the Manager Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication December 31, 
1980. This paper part the Journal the Irrigation and Drainage Division, Proceedings 
the American Society Civil Engineers, Vol. 107, No. IR3, September, 
1981. ISSN 0044-7978 /81 /0003-0247 
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and planning consistent with projected climate changes anticipated during the 
next years. 

Projected Climate Change.—The increasing use fossil fuels probably 
the most important economic and environmental issue facing mankind. During 
June 19-20, 1980, workshop was held the Institute for Energy Analysis, 
Oak Ridge, Tennessee, Change and Food Production Nonindus- 
trialized The workshop was sponsored the American Association 
for Advancement Science, AAAS, part project jointly sponsored 
AAAS and the U.S. Department Energy, DOE, and 
Societal Consequences Possible Carbon-Dioxide Induced Climate 
Some participants from eight countries developed draft research proposal 
for submission DOE AAAS. 

Although many scientists and meteorologists have not yet accepted the 
conclusion that temperatures will rise significantly during the next years, 
the workshop did not debate question this conclusion but rather attempted 
develop proposal for research that would maximize possible benefits 
agriculture and food production while also minimizing the undesirable effects. 

indicated DOE (7), Manabe and Wetherald (6), and Wigley, Jones 
and Kelly (9) the use fossil fuel increases the current rate, CO, 
the atmosphere will approximately double within years. Wigley, Jones, and 
Kelly (9) estimate that doubling CO, the atmosphere will produce expected 
increase global mean temperature about They show that associated 
decrease rainfall will occur for most North America with the change 
warmer years. Manabe and Wetherald (6) developed simple general circulation 
model with seasonal change insolation and simplified interaction between 
cloud and radiative transfer. The computational domain divided into two 
equal areas continent and ocean. comprises 120° longitude and 90° 
latitude. their model, predicted precipitation over the continent decreases 
between latitude 37° and 45° with doubling Evaporation and evaporative 
demand increases nearly all latitudes. Precipitation minus evaporation has 
the largest negative values between 37° and 47° latitude. these latitudes, 
the mean projected temperature increase for 100% increase atmospheric 

The model, because its many simplifications and idealizations, should 
used with caution. Quantitative details may have significant error either 
side more realistic values. was developed provide preliminary insight 
into the possible climatic response the increase CO, concentrations and 
establish foundation for future global model having realistic geography. 

The Department Energy (7) cites 1975 publication Manabe and Wetherald 
(5) and indicates increase mean temperature 2.8°C 37° 
47° latitude. They show most the United States becoming drier. The 
following quoted from the Department Energy (7): 


the sense the scientific community that carbon dioxide from 
the unrestrained combustion fossil fuels potentially the most important 
environmental issue facing mankind. Current predictions call for doubling 
atmospheric carbon dioxide early 2030. Climate models using 
these elevated levels, predict the possibility significant dislocations 
the global distribution climate. Should such perturbations the distribu- 
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tion global temperature, cloudiness, precipitation, and wind occur, then 
clear that there could major changes the global distribution 
agriculture the extent and composition ecological systems, 
the biological and physical characteristics the oceans, and the extent 
the cryosphere. The costs benefits such extensive environmental 
changes balanced international societal system are difficult estimate 
but could historically unprecedented. therefore the utmost 
importance that the technical community resolve the scientific uncertainties 
concerning the effect carbon dioxide climate and the cascade 
environmental and societal impacts resulting from change climate 
quickly possible. 


Temperature Change.—The estimation future mean temperatures the 
earth’s surface requires assumptions the amount forest clearing, fossil 
fuel production and use, and the amount carbon dioxide that will remain 
the atmosphere. The possible margin error these estimates considerable. 
The model prepared Manabe and Wetherald (6) indicates rather modest changes 
2°C less near the equator and maximum increase temperature 
about 83° latitude. For doubling CO, the maximum increase and 
for quadrupling CO, shown 15° 

The estimated increase not completely linear but for the range latitudes 
65° the temperature increase, °C, can estimated reasonably 
well from latitude, LA, using linear equations. For doubling carbon dioxide 
the equation is: 


and for quadrupling carbon dioxide 


From these predicted temperature changes seems probable that doubling 
carbon dioxide the atmosphere would shift climatic zones for crop adaptation 
about 250-300 (400-480 km) further north. The frost-free growing season 
would increased about one month for latitudes exceeding 30°. 

increase CO, will also increase photosynthesis. Higher temperatures 
would increase photosynthesis for some crops, particularly the crops 
those crops with four carbon atoms the carbohydrate rings. crops include 
corn, sugar cane, sorghum, millet, and some the grasses and forage plants. 
Water requirements would increased due the higher temperatures and 
the longer growing seasons. 

Estimating Potential Evapotranspiration, purposes this study, 
potential evapotranspiration, ETP, taken index water requirements. 
The question arises how ETP will vary with increasing temperature. ETP 
also influenced radiation, and, latitudes 37° 47°, rainfall projected 
decrease. The decrease not great, but, since the amount radiation reaching 
the earth depends upon extraterrestrial radiation and upon how much absorbed 
and scattered while passing through the atmosphere, radiation may significantly 
increase. Absorption and scattering radiation largely produced clouds 
and water vapor. With decrease rainfall and cloudiness, increase 
incident solar radiation can anticipated. However, this increase should not 
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large, and, since would quite difficult estimate the magnitude 
the increase, not taken into account this study. 

The California Department Water Resources (8) selected Class pan 
evaporation from pan located large, irrigated pasture the most practical 
and reliable means estimating ETP. This decision was reached based upon 
several years data collection from lysimeters located irrigated pastures 
and comparisons with methods estimating ETP from climatic data. The 
Technical Committee Irrigation Requirements the Irrigation and Drainage 
Division, ASCE (1), compared methods for estimating ETP using world-wide 
lysimeter measurements. The Committee report states careful standard- 
ization pan environment maintained, and where strong, dry-wind conditions 
occur only occasionally, the use pans may reliable any the 
methods discussed However, the data from California indicate 
that the use large, irrigated pasture eliminates the need for the qualifications 
stated the ASCE Committee. 

Using five-day sums grass lysimeter evapotranspiration, ETL, temperature 
degrees Fahrenheit, incident solar radiation, RS, relative humidity, 
RH, and wind movement, from eight years daily data collected Davis, 
California, was found that predicted 94% the variance ETL, 
0.94, for five-day sums. percent predicted 9.2% the remaining 
variance and wind only 1.5%, showing that adding corrections for wind 
and relative humidity equation could predict nearly 95% the variance 
ETL. Since other factors did not significantly improve the prediction, 
equation was developed using and only. The equation can written: 


which ETP and are the same units equivalent water evaporation. 

Eq. has been compared with lysimeter evapotranspiration and pan evaporation 
for eight locations representing wide variety climatic conditions California 
using data from the California Department Water Resources (8). There 
little significant difference between accuracy predicted values from well 
standardized Class pan located large irrigated pasture and those from 
Eq. means for estimating ETP. 

Although the ASCE Committee has refused make the lysimeter and other 
data available for others use check methodology, data obtained independently 
from four locations indicate that Eq. superior means estimating 
ETP any the procedures tested the Committee. Hargreaves (2) 
compares four locations with results from Kohler, Penman, Christiansen, Jensen- 
Haise and Stephens-Stewart. 

Hill (4), review the methods tested the ASCE Committee, describes 
all the methods being somewhat site specific. crop coefficients are 
generally applicable, the method for estimating ETP should produce values 
such that ETP will zero when ETL zero. The relationship between ETL 
and ETP should approximately linear with one one slope. Eq. seems 
superior the methods tested this criteria used for evaluation. 

Most the equations evaluated had variable intercepts and variable slope 
relationships when applied differing lysimeters, indicating that either the 
equation crop coefficients require local calibration. Eq. not site specific. 

Mean temperature used most evapotranspiration estimating methods. 
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ETP principally daylight phenomenon. high elevations, the difference 
between maximum and minimum temperatures much greater than low 
elevations. This introduces bias the use most equations, including Eq. 
when used for high elevation locations. Also advective energy often accounts 
for significantly larger percentage ETL higher elevations. Some correction 
Eq. may desirable for higher elevations particularly the conditions 
favor the advective transfer energy. 

major limitation the evaluation estimating methods the paucity 
reliable solar radiation measurements. Radiometers vary type and accuracy. 
Some are not adequately calibrated checked for drift. humid areas, rain 
and condensation sometimes reduce radiation readings significantly below actual. 
Preoccupation with including various weather elements predictive equations 
will produce limited benefits long one the dominant values temperature 
and radiation frequently measured with accuracy plus minus 10% 
and sometimes 20% more. 

Hargreaves (3) used the climatic mean values temperature and estimated 
radiation from data published the World Meteorological Organization, WMO, 
estimate ETP for 644 world-wide locations. was estimated from extraterres- 
trial radiation, RA, and percentage possible Sunshine, The equation can 
written: 


When values were not available they were estimated from relative humidity, 
RH, using the equation: 


which constant that requires local calibration for each climatic zone 
and for the hours and methods measuring relative humidity. 

Water Use Changes.—For Mexico, the United States, and Canada, the estimated 
values ETP were graphed function latitude, LA. For the latitudes 
15° 65° North, annual ETP for low elevation can 
approximated the equation: 


From Eq. and the estimated temperature increases resulting from the doubling 
the concentration carbon dioxide, estimate was made the increase 
ETP due higher temperatures. The results are presented follows: 


Latitude degrees 
percentage 


For latitudes 30° 60°, there added increase ETP due the lengthening 
the growing season. The projected temperature change would shift present 
mean temperatures for latitudes 30° 60° about two degrees further from 
the equator and increase the length the frost free season about days, 
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resulting the following percentage increases the length the growing 
season: 


Latitude degrees 
Increase growing sea- 


assumed that the percentage increases ETP for increasing temperature 
and the percentage increases length growing season can added indicate 
total resulting increase ETP, then the total increases are follows: 


percentage 


Problem Area.—For latitudes 37° 47° Central United States, the Manabe- 
Wetherald model indicates decrease precipitation and therefore decrease 
soil moisture available for plant growth. Considering both the higher temperature 
and the increased length growing season seems probable that about 
percent more water would required during the growing season order 
fully utilize the increased potential for production. 

42° north latitude, average annual potential evapotranspiration low 
elevations about 1,000 mm. 30% increase equivalent 300 
depth 11.8 in. These estimated changes are significant since the major portion 
the present corn and livestock production the United States located 
the potential problem area. Both domestic and world food requirements dictate 
need minimize losses production from this rich agricultural region. 
would seem that immediate attention needs given studies designed 
make this potentially large water deficit through large scale mass water 
transfer over long distances. 

The climate model used and the ETP methods presented herein are not precise 
and the climatic changes may significantly greater considerably less than 
predicted. AAAS-DOE workshop, was generally accepted that Central 
America would benefit from increased rainfall the future. However, the amount 
rainfall Central America presents less problem than temporal and 
physical distribution. For much the Central American Isthmus, the rainy 
season six months, contains period July, August, 
both, referred the little dry season the 
regional meeting hydrologists and meteorologists representing the six 
countries held Tegucigalpa, Honduras, July, 1980 (regional meeting 
Comite Regional Recursos Hidraulicos, Apartado 718 Tegucigalpa, Honduras), 
attempt was made evaluate the effect CO, induced temperature 
increase agricultural production. Based conclusions the most complete 
analysis, increase mean temperature would increase annual rainfall but 
would also increase the amounts during periods that are now excessive and 
would extend and intensify the little dry season. These changes would increase 


the need for all kinds water control, including irrigation, drainage, and flood 
control. 
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Benerits 


increase CO, will increase potential plant growth. Higher temperatures 
will make larger yields possible for some crops. Good lands northern latitudes 
will become more productive. water can made more available and used 
more efficiently, the lands receiving less rainfall will have much greater potential 
for agricultural production. However, major problem may the level 
human and animal tolerance heat. 

the latutides 37°-47°, the projected temperature increase from quadrupling 
carbon dioxide concentrations 6°-7° average increase about 12° 
order for food production potential fully benefit from this change, water 
requirements would greatly increased over present levels. Plant and animal 
genetic material would need developed modified incorporate adjust- 
ments the higher temperatures. 

Another possible benefit potential climatic change studies that clearer 
understanding the probable increases temperature will speed the 
technological advances required for developing alternatives the use fossil 
fuels. Although potential benefits climatic change seem large, there also 
the question how much change can reasonably tolerated and what 
point negative aspects outweigh potential benefits. 


The AAAS-DOE sponsored workshop selected nine important research topics 
with various subheadings. These were listed order priority. The following 


topics relating irrigation and water management are listed follows with 
their original priority numbers: 


The management water rainfed agriculture and irrigation frequency. 
Better understanding climate/production relationships and collection 
and use necessary information. 

Direct effects CO, photosynthesis, water use, and nitrogen fixation. 
Soil temperature and moisture relationship and assessment. 


Variability agricultural production from year year and from one location 
another influenced more water than any other factor consideration. 


For this reason, the working group assigned first priority research related 
water management. 


Summary 


The concentration carbon dioxide the atmosphere may double within 
the next years result increased use fossil fuels. Temperatures 
are projected rise. Photosynthesis and plant growth will generally increased. 
Between latitudes 37° 47°, rainfall indicated decreasing. Potential 
evapotranspiration these latitudes may increase much 30%. Future 
agricultural production will become increasingly dependent upon the adequacy 
water supplies for irrigation. Climatic zones suitable for specific crops may 
shift 250-300 further from the equator. 
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With increasing world population and major changes location and amount 
future irrigation land with potential for agricultural production, will 
necessary revise priorities for water resource development and reassess 
the role that the United States can and should have major producer and 
exporter food. 

Toleration temperatures associated with doubling carbon dioxide con- 
centrations may not present unmanageable problems for mankind for plants 
and animals. However, many adjustments will necessary. Contemplation 
higher temperatures associated with quadrupling carbon dioxide concentra- 
tions raises questions the tolerable upper limit for plants, animals and 
humans. This requires much further research and analysis. 

some levels change, appears that modifications induced increased 
carbon dioxide concentrations can produce fairly substantial benefits, providing 
water made available deficient areas. Since major man-induced changes 
are projected within the lifetimes individuals now living, there real need 
act quickly document present climatic conditions more carefully, reduce 
the margin error estimating the changes and take the necessary steps 
maximize benefits and minimize the possible dislocations and undesirable 
effects. 

order model the effects climate change world food production, 
improved information the water balance required. Rainfall needs 
estimated during periods that correspond the capacity soils store useful 
water for crop growth and potential evapotranspiration. Reliable estimates 
ETP require accurate measurements estimates solar radiation. Therefore, 
highest priority needs assigned developing world-wide information 
present conditions precipitation, solar radiation, and water balance, all 
which provide good spacial coverage. Improved modeling future water balance 
conditions can then undertaken. Without such improved information and 
models, satisfactory estimates the probable impacts climate change 
world food production will impossible. 
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The following symbols are used this paper: 


ETL grass lysimeter evapotranspiration; 
ETP potential evapotranspiration, usually mm; 
constant used estimating sunshine percentage from relative 
humidity; 
latitude degrees (in this report North latitude); 
extraterrestrial radiation equivalent water evaporation; 
relative humidity percent; 
incident solar radiation equivalent water evaporation; 
percentage possible sunshine; 
mean temperature degrees Fahrenheit; 
temperature increase (resulting from increase CO, concentra- 
tions the atmosphere); and 
wind movement distance per unit time. 
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ADVANCE RATES UNDER SURGE 
FLow SYSTEMS 


Niel and Gregory Poole* 


INTRODUCTION 


Considerable attention has been given the advance phase surface irrigation 
because its bearing intake opportunity times between the two ends 
the irrigation run. Surface irrigation design and management objectives are 
generally complete the advance phase the irrigation quickly possible 
minimize these differences intake opportunity time. accomplish this, 
large streams have been recommended for the advance phase with cutback 
streams runoff recovery systems minimize the runoff and deep percolation 
losses during the intake phase. The ideal system would capable automation, 
provide instantaneous advance, and prevent runoff during the intake phase. 
Unfortunately, most surface irrigation systems, particularly sloped furrow sys- 
tems, not have the physical resources needed for achieving these conditions 
without substantial capital investment. 

Surge flow for automatic furrow irrigation was first suggested Stringham 
and Keller (6) improved method automating cutback furrow irrigation. 
Since most simple automatic irrigation valves can only turn water and off, 
they concluded that would simpler cycle valves reduce the flow 
rate time basis than partially close the valves achieve the cutback 
stream. cycling flow system complete with valves and automatic controller 
was developed and successfully tested accomplish pulsed application water 
the furrows. The field testing the flow equipment indicated that cutback 
could achieved almost any desired average stream size using time cycles 
from few seconds hours duration. Perhaps more importantly, was 
found that the cycled flows had major impact furrow advance rates 
well the temporal and spatial intake variability within the field. 
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Surce 


The intermittent application irrigation surface irrigation systems 
distinctly different than with either sprinkle trickle systems. Karmeli and 
Peri (4) described the pulsing sprinklers achieve reduced time averaged 
application rates. furrow system, however, pulse water introduced 
the furrow inlet subject the various resistance and inertial forces 
moves along the furrow, and therefore elongates and slows down. important 
result these dynamic factors that the advancing tip moves more rapidly 
than the receding edge. Thus, when the intervals between water inputs shorten, 
the pulses eventually merge together some distance along the furrow. fact, 
during the early testing Stringham and Keller (6), the flow the end 
the furrows was relatively steady. The behavior pulsed furrow inflows 
better characterized surge hydraulics, defined Chow (3). The term 
surge flow was adopted and still conceptually correct when time 
intervals between surges great enough that individual impulses not 
interact. 

Surge flow creates series and off conditions constant variable 
time spans the furrow inlet. The cycle time has been defined the period 
required for complete on/off cycle, i.e., the time between the beginning 
one surge the beginning the next. The cycle time may any desired 
duration and can vary from few seconds hours, but present experience 
indicates cycle time min-60 min. The cycle ratio defined 
the ratio the time cycle time, e.g., with conventional irrigation, 
furrow set might considered have 6-h cycle time with 
off time, making the cycle ratio equal one (continuous flow). This same 
irrigation might accomplished with the same furrow stream and irrigation 
time using 30-min cycle time and cycle ratio one-third, 
min and min off each cycle. The applied water would reduced 
one third, reducing the time for each furrow from but allowing 
more convenient total irrigation time. The number furrows per 
set would increased decreased proportionately. 


Field studies were carried out during the 1979 irrigation season test the 
characteristics surge flow irrigation furrows with particular emphasis 
the advance phase. The third writer (1) used instantaneous furrow streams 
gpm (0.63 L/s), gpm (1.26 L/s), and gpm (1.89 L/s) with cycle 
ratios one for the gpm (0.63 L/s) inflow, one half for the gpm (1.26 
L/s) stream and one third for the gpm (1.89 L/s) stream, resulting 
equal quantity water being applied each furrow over given time period. 
The cycle time for these tests was min. The fourth writer (5) used surge 
flow furrow streams gpm (1.26 L/s) with cycle times min, min, 
min, and min. each case the cycle ratio was one half, making the 
time averaged flow stream gpm (0.63 L/s) and equal continuous flow 
gpm (0.63 L/s) furrows irrigated the same time for comparison. 

The studies were carried out the Utah State University Dairy Research 
Farm located about mile north the Utah State University campus. The 
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soil classified silt loam and the site was planted corn with 
in. (0.76 row spacing. The furrows were 600 (183 long and had 
average slope the direction irrigation 1.46%. Planting and cultivation 
were planned that tractor wheels always ran along the same furrows, resulting 
compactive effort every other furrow. The compacted furrows were 
designated ‘‘wheel and those not compacted the tractor wheels 
were designated The tools used form the furrows 
before the first irrigation covered the wheel tracks and left all furrows with 
similar roughness and cross section. Prior making test run, 100 (30.5 
stations were established along the furrows aid measuring the advance 
rate. addition, soil moisture samples were taken from depths 0.5 ft, 1.0 
ft, 2.0 ft, and 3.0 (150 mm, 300 mm, 600 mm, and 900 mm) for every other 
test furrow. 

Using the equipment described Stringham and Keller (6), water was applied 
randomly selected furrows the cycle times and cycle ratios outlined 
previously. During each test, the time for water advance each station 
along the furrow was recorded. Furrow runoff was also measured that intake 
could determined from the inflow-outflow data. Since very little information 
exists regarding surge flow, several exploratory runs were made before choosing 
the cycle times and cycle ratios decided upon. 

Several replications each treatment were conducted over the field and 
was initially assumed that spatial variability could neglected. Consequently, 
all the advance data were first fitted power advance function and then 
treatment produce single relation each cycle time and 
cycle ratio combination. The mean relation was determined numerically 
averaging the logrithmically transformed power advance coefficients. 


The effects surge flow irrigation were most apparent during the first irrigation 
when the hydraulic and infiltration characteristics were extreme. Fig. shows 
the advance data for the first irrigation non-wheel furrows July 1979 
(1). The average flow rate gpm (0.63 L/s) for continuous flow treatment 
and three surge flow treatments utilizing min cycle time and cycle ratios 
shown. can seen that the continuous flow treatment required almost 
order magnitude more time complete the advance than the surge flow 
irrigated furrows. Although the difference between the continuous and surge 
flow treatments statistically significant the confidence level, the 
differences among the surge flow treatments were not. Fig. shows the results 
the same tests for the second irrigation. The advantages surge flow were 
substantially reduced. The differences surge flow advance rates are not greatly 
different but the continuous flow produced much faster advance. 

each case previously noted for the third writer (1) data, the time average 
furrow inflows were equal gpm (0.63 L/s). eliminate the effects 
variable instantaneous flow rates, the fourth writer (5) introduced gpm (0.63 
L/s) into the continuous flow furrows and into each surge flow furrow but 
varied the cycle time. The cycle ratio was fixed one half. Again, the most 
pronounced effects surging the furrow inflows occurred during the first 
irrigation the wheel furrows (Fig. 3). can seen that the advance under 
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FIG. Rates under Continuous and Variable Cycle Ratio Surge Flow 
Treatments Newly Formed Non-Wheel Furrows 0.305 


FIG. 2.—Advance Rates for Continuous and Variable Cycle Ratio Surge Flow Treat- 
ments Non-Wheel Furrows During Second Irrigation 0.305 
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FIG. 3.—Advance Rates Newly Formed Wheel Furrows under Continuous and 


Two Variable Cycle Time Surge Flow Treatments 0.305 


C - Continuous 1wFLOW 

> - S-MINUTE CYCLE INFLOW 
10-MINUTE CYCLE INFLOW 
20-mMINUTE CYCLE INFLOW 


ADVANCE TIME (MIMUTES) 


DISTANCE 
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FIG. 5.—Advance Relations for Various Continuous Flow Furrow Irrigations 1979 


FIG. Relations for Surge Flow Furrow Irrigations 1979 0.305 
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surge flow conditions was approximately four times faster. Results from the 
first irrigation the non-wheel furrows are shown Fig. After the first 
irrigation when infiltration differences should less noticeable the field 
when wheel compaction reduces these differences mechanically, the advance 
under surged conditions were much closer the rates under continuous flow. 
appears that introducing furrow flow and then letting recede and infiltrate 
before the next surge passes significantly affects the permeability the furrow 
perimeter. Given the tendency for newly formed furrows have widely different 
intake characteristics than previously irrigated furrows, also appears that surge 
flow systems may significantly reduce these differences. 

These effects surge flow practices infiltration rates are probably the 
most important aspect this new irrigation methodology. illustrate these 
effects further, the data the fourth writer (5) were grouped continuous 
treatment Fig. and for the min cycle Fig. While the shortening 
the advance time evident, previously noted, more interesting observation 
the total variation between the two treatments. Over the season, over the 
field, and among replications, advance times ranged from 270 min-3,490 min 
for the continuous flow furrows, irrespective compactive history. For the 
min cycle, the range varied from min. other words, the variation 
under continuous flow conditions was approximately factor 13, whereas 
for the min cycle was only factor about two. further testing 
the surge flow concept verifies significant reduction the temporal and spatial 
variability intake, surge flow systems will amenable much more precise 
design and operation strategies than for conventional furrow irrigation. 


The research reported herein concerns single field condition and single 
time average furrow inflow rate. Instantaneous furrow inflows were varied 
conducting surge flow tests having variable cycle ratios. all trials date, 
the advance phase was accomplished with less water using surge flows than 
with continuous flows. The effects surge flow were most pronounced during 
the first water application newly formed furrows and were also more 
pronounced non-wheel furrows. comparison data aggregated treatment 
but not time location within the field showed that differences advance 
rates under surge flow conditions were substantially less than under continuous 
flow conditions. 

conducting these tests and evaluating the data, it-has become obvious 
that surge flow alters the basic intake characteristics the furrow. theorized 
that cycled movement water over the soil surface accelerates the development 
thin surface seals the bottom the furrows. Apparently, during the interval 
between surges, the bottom soil consisting dispersed fine materials and 
lubricated water are compacted tension forces which build the 
soil the furrow continues drain. 

The limited research date indicates that surge flow has the potential 
powerful management technique improve the application uniformity and 
efficiency furrow irrigation. The volume water necessary complete 
the advance phase can substantially reduced using surge flow over the 
conventional practices. Likewise, the intake opportunity time between the head 
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end and the far end the run can equalized, resulting more uniform 
water distribution. Finally, the possibility programming changes the cycle 
ratio achieve cutback can effectively eliminate runoff with its water loss 
and irrigation return flow problems. 

evident, however, that much additional study required establish 
the design criteria for the wide range soils, length run, size stream, 
and other conditions surface irrigation. The effect surge flow intake, 
advance, uniformity, and efficiency must determined. Theoretical work 
the dynamic geometry the surge and its relationship intake and advance 
necessary. 


ACKNOWLEDGMENT 


The study reported herein part research project No. 794, Utah Agricultural 
Experiment Station, Logan, Utah. 


Allen, L., Rates Furrow Irrigation for Cycled thesis presented 
the Utah State University, Logan, Utah, 1980, partial fulfillment the 
requirements for the degree Master Science. 

Bishop, A., Surge Utah Science, Vol. 41, No. 1980, pp. 
60-64. 

Chow, T., Open Channel Hydraulics, McGraw-Hill Book Co., Inc., New York, 
N.Y., 1959. 

Karmeli, and Peri, G., Principles Pulse Journal the Irrigation 
and Drainage Division, ASCE, Vol. 100, No. IR3, Proc. Paper 10831, Sept., 1974, 
pp. 309-319. 

Poole, G., and Advance under Surge Flows Furrow Irrigation,’’ thesis 
presented the Utah State University, Logan, Utah, 1980, partial fulfillment 
the requirements for the degree Master Science. 

Stringham, and Keller, J., Flow for Automatic Irrigation,’’ presented 
the July, 1979, ASCE Irrigation and Drainage Division Specialty Conference, held 
Albuquerque, N.M., pp. 132-142. 


SEPTEMBER 1981 


DEVELOPMENT SOLUTIONS 
FOR LEVEL-BASIN DESIGN 


and Allen 


INTRODUCTION 


Level-basin irrigation has grown rapidly popularity over the last decade. 
There are several reasons for its popularity: (1) Uniform application water 
and associated high application efficiencies, which reduce farm water require- 
ments and decrease nutrient and fertilizer loss from excessive deep percolation; 
(2) elimination runoff, which reduces soil erosion and eliminates water reuse 
systems; and (3) better control the irrigation water, thus reducing labor 
requirements. The development laser-controlled land leveling equipment has 
greatly improved the achievable efficiencies for level basins. 

level-basin irrigation, water enters the field one end corner 
and spreads over the basin. The driving force the water surface slope—the 
higher the flow rate, the greater the flow depth and water surface slope, and 
the faster the rate advance. Once the entire basin inundated the water 
surface levels out. The irrigation over when all the water infiltrates into 
the ground (recession). This recession occurs about the same time over the 
entire basin. The difference opportunity time from one part the field 
the other essentially the difference advance time. Since the rate 
advance increases flow rate increases, the uniformity application will 
also increase with flow rate. 

Level-basin irrigation well suited low intake soils since the difference 
infiltrated depth for given difference opportunity time relatively small. 
Also, advance across the basin faster for low intake soils than for higher 
intake soils since more water remains the surface, increasing the water surface 
slope. This reduces the problem balancing stream size and application time 
commonly associated with sloping-border irrigation. This paper addresses proce- 
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dures for obtaining reasonable estimates distribution uniformities for wide 
variety flow rates, lengths run, infiltration characteristics, and flow resistance 
for the design and management level basins. 


Basin 


designing level basins, there are many options consider depending 
imposed limitations stream size, basin size, soil infiltration characteristics, 
and crop management considerations. many cases, there maximum practical 
stream size depending the delivery source and field turnouts. Erosion becomes 
problem with high flow rates the water must spread out over large 
area using several small turnouts, erosion control structures, both. Maximum 
water depth can also impose limitation stream size, particularly with narrow 
basins row crops. Basin size often limited physical constraints, such 
roads, property lines, canals, etc. 

The Soil Conservation Service (SCS) has developed design manual for border 
irrigation (20). series design charts are presented for level-basins and sloping 
borders. The charts presented for level basins relate unit flow rate and length 
run for various application efficiencies. range Manning roughnesses, 
soil infiltration families, and net depth application are presented the charts. 
Kruse defined the term application efficiency the ratio the average 
depth the irrigation water stored the root zone the average depth 
irrigation water applied (13). the irrigation completely fills the root zone, 
the application efficiency for unit width border 


which the field length; the flow rate per unit width; the 
application time time cutoff; and the soil moisture deficit the 
depth required fill the root zone. also useful define the distribution 
uniformity for blocked end borders basins 


which the minimum net depth infiltrated. The relationships between 
Z4, DU, and are shown Fig. where the gross average 
depth applied, When the minimum depth infiltrated less than 
the desired depth (z, z,), then [Fig. and the entire field 
has been over irrigated. When z,, then and maximum 
for adequate irrigation [Fig. However, when z,, then 
and part the field has been underirrigated. Note that for this case, 
not easily calculated [Fig. similar analysis can made for irrigation 
efficiency, which includes other beneficial uses water. 

The SCS solution assumes that and the soil moisture 
deficit just met the end the field. The solutions developed this paper 
use DU, since more straightforward for all three cases. Conceptually, the 
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case Fig. the most desirable, and for design purposes can determine 
and use the potential general, the distribution uniformity 
dependent the irrigation system design. Application efficiency depends 
both design and irrigation management and individual irrigation and 
seasonal parameter. The distribution uniformity can thought the upper 
limit application efficiency (when not considering deficit irrigation net 
soil moisture depletion over the growing season). 

well designed systems, both sloping (without pumpback) and level, higher 
efficiencies are usually attained the incoming stream cut off before the 
entire field has been covered with water. However, the SCS solution method 
for level basin design assumes that advance complete before cutoff. They 
also make several other arbitrary assumptions regarding the friction slope and 
water depths (5). Also, the SCS intake families, principle, couple the magnitude 


DISTANCE 


DEPTH 


DU< 


FIG. 1.—Profiles Infiltrated Water: (a) For Overirrigation Throughout Basin; (b) 
For Adequate Irrigation with Excess Overirrigation; and (c) For Slight Underirriga- 
tion over Part Field 


and shape the infiltration function which not always represent the true 
infiltration the soil. 

Merriam (14), well others, emphasized that the time required infiltrate 
the desired depth application soil moisture deficit key factor designing 
border-irrigation systems. For analyzing distribution uniformity, more useful 
use the minimum opportunity time, which corresponds z,,. 

Once the basins have been designed and built, management the system 
consists determining the application time necessary apply desired volume 
water (or gross depth application). The application time time cutoff 
can found from provided that the flow rate known and accurate 
estimate the distribution uniformity obtainable. 
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The appropriate factors which govern the flow irrigation borders basins 
are the inflow rate per unit width the flow resistance (Manning n), time 
cutoff application time field length (L), field slope (S,), the field 
end conditions (blocked end open ended), and the infiltration characteristics 
the soil. 

Field evaluations all the possible combinations the aforementioned factors 
would virtually impossible. For this reason, desirable use mathematical 
model the system determine the results different combinations these 
variables. Many models have been proposed and developed describe the 
flow water irrigation borders. Some are based empirical approaches, 
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DIMENSIONLESS INFILTRATION TIME, 


FIG. Distribution Uniformities Representing Variations Infiltration 
and Soil Moisture Deficit, 0.5 


volume balance techniques, simplified hydraulic models 
These methods depend normal depth the surface stream, thus they 
not apply the level basin case. Other models are based exact mathematical 
equations for flow open channels, i.e., the Saint Venant equations (1,10,12). 
While highly accurate, these models are time consuming and expensive run. 
Recently, new model, based the Saint Venant equations 
has been developed (19). This model accurate and relatively 
reliable, but inexpensive run. Verification the model was presented 
Strelkoff and Katopodes (10) and Clemmens (3,4). 
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Given the unit inflow rate and cut-off time, and the length, roughness, and 
infiltration characteristics the basin, the mathematical model capable 
predicting the advance and recession the surface stream, and thus, the 
post-irrigation profile infiltrated water over the length the basin. These 
advance and recession data were used develop the uniformity curves presented 
this paper. Given the minimum depth application z,, the distribution 
uniformity, DU, can computed (as can other measures the merit 
given choice input variables): which and 
are the constant and exponent power infiltration function which will 
defined later section. The variable not independent variable; 


DIMENSIONLESS FLOW RATE, 


FIG. 3.—Level-Basin Distribution Uniformities Representing Variations Flow Rate 
and Field Length, 0.5 


computer runs, with systematic variation each the independent variables 
(input parameters) would show the effect distribution uniformity. However, 
astronomical number runs and graphs results would required. The 
number runs and output graphs reduced manageable size using 
dimensionless notation. will shown the following sections, that terms 
depth infiltration relative characteristic depth; the time required 
infiltrate depth z,, relative characteristic time. This function, plotted 
Fig. for 0.5, can used show how the minimum infiltrated depth 
and infiltration characteristics the soil (as well field length, roughness, 
and inflow rate) affect the distribution uniformity. 
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Alternately, will shown that L*, which shows, with 
different set dimensionless terms, the effect inflow rate and field length 
(as well minimum depth, field roughness, and infiltration characteristics) 
the distribution uniformity. This function plotted, for 0.5, Fig. 
Illustrative examples show how these curves can used determine suitable 
values inflow rate, cut-off time, and basin dimensions. 


Basic 


The second and first writers (18) presented detailed dimensional analysis 
the variables associated with both sloping-border and level-basin irrigation. 
The dimensionless, differential form the Saint Venant Equations for level 
basins, resulting from this analysis, 


the dimensional form these equations, flow rate per unit width; 
distance down the border; surface water depth; infiltrated depth; 
time; infiltration opportunity time; and relative resistance flow, 
drag per unit plan area surface stream per unit weight water (friction 
slope times y). Starred values for (and and are dimensionless 
values these variables that are defined terms reference variables, 
three variables are the result the nondimensionalization, namely 


Boundary conditions for the solution Eqs. and for basin irrigation are 
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which the unit flow rate the head the basin; 
cutoff application time; the advance distance; and 
length. 

The solution Eqs. and depends upon values chosen for K*, 
D*, and L*. The solutions for level-basin advance presented the 
loss generality. Furthermore, was postulated that Finally, the 
characteristic infiltration depth was defined z(T). The solution was 
governed and additional parameter resulting from selection the 
infiltration function. The number parameters that can arbitrarily set and 


the number parameters governing the solution were considered detail 
the second and first writers (18). 


the field 


Systems 


The solution presented here assumes power infiltration function the form 


and form the Manning equation for roughness 


7/3 


which units coefficient (1.0 the metric system and 1.486 
the English system). Similarly, Use Eq. adds 
one additional governing parameter, which must specified the dimension- 
less solution (18). 

The dimensionless equation for distribution uniformity 


with and the only other possible governing parameters governing DU. 
The solution Eqs. and 8-12, with K*, D*, and 
all fixed, leads distribution infiltrated depths resulting from the advance 
and recession the irrigation stream. Thus, for each set input conditions, 
value for can obtained and calculated from Eq. 15. The resulting 
solution for would governed and two other dimensionless parameters. 
For the dimensionless system, used the first and second writers (6), which 
values for K*, and L*. Alternately, the solution for could presented 
function z*, and first glance, appears that this direct solution 
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From Eq. 15, with 


Obtaining the dimensional equivalent this form dimensionless solution 
would require that and all known priori. This impractical 
for design and management purposes, since not known advance. 
would more practical were dependent, rather than independent, 
variable. Thus, for useful solution, should not one the fixed parameters. 

The possible choices for alternate fixed parameter replace 
result the solution for advance and recession and not known 
priori. Thus, trial and error solution would necessary 
Selecting was also considered the first and second writers (6). 
They did not seiect it, since advance after cutoff could displaced more 
desirable, selection logical choice. avoid confusion, the 
following subscript notation was adapted distinguish between dimensionless 
systems: the work the first and second writers (6), and 
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can replaced the list independent dimensionless variables affecting 
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the new independent variable since This replacement 
has two advantages. and (or and have more physical significance 
than and (or and relation surface irrigation. Further, 
replacing the characteristic depth infiltration Y,, the parameter 
(and Eq. 22) removed from the calculation reference variables. Thus, 
the solutions for with different values can compared without 
recalculating the reference variables and sensitivity analysis can made 
direct observation the solution curves. 

Advance and recession curves for this dimensionless representation could 
Then, the solution for could displayed terms and 
Computationally, there one drawback this dimensionless solution. The 
solution technique for solving the zero-inertia equations uses discrete time steps 
for the computation distance increments during advance. Advance considered 
complete when L*. Obtaining solution that results value 
exactly one unlikely. Thus, would more precise compute the solution 
with and transform the results, considered follows. 

For given dimensional set values for and z,, there 
one and only one value DU. Thus, once the value for determined 
from the solution the governing equations, can displayed any form 
long the same for that same set dimensional values. Since 
the relationships between the reference variables must maintained, Eqs. 5-7 
and must also still hold. Thus, indirect solution can obtained. The 
possible choices available for reference variables (two which 
can fixed) are and Two these combination have been 
the remaining combinations, only z,, practical significance. 
This dimensionless system has been denoted the subscript Specifying 
also gives since z(T). Thus, the dimensionless system 
variables from 5-7 gives 


and from 


The solution can presented terms and Also, with 
the previous (B) dimensionless system, not involved the calculation 
the reference variables. 
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The governing equations for flow rectangular level basin (Eqs. 
and 8-12) were solved for the postirrigation profile infiltrated water with 
ranging from (depending The distribution uniformity was computed 
with Eq. for each K*, combination and the resulting value 
set values pertaining given values DU, i.e., 95, ... 50, was 
then determined for each value with interpolating polynomials. These 
results were then transformed represent solutions terms the two 
The transformation was made point point using Eqs. 16-26 
with all dimensional variables held constant. simplify the computations, 
equations, were developed relating dimensionless values one system those 
another. Curves through these points were then smoothed yield the curves 
presented Figs. and Fig. can used determine the effect the 
infiltration function the distribution uniformity. The axes Fig. are made 
dimensional multiplying values along the axes and from 
and (similarly for Fig. with and from and 25). Fig. 
can used determine the effect unit flow rate and basin length 
distribution uniformity. Fig. resembles the level-basin design curves presented 
the SCS their border irrigation handbook (21). Unlike the SCS graphs, 
these curves allow for continuous, rather than discrete, variations Manning 
roughness, net infiltration depth, and the corresponding required opportunity 
time. Both dimensionless representations can used for evaluating existing 
systems, for designing new level-basin systems, and for developing management 
aids for actual operation. These curves represent solution for value 
0.5 only. complete solution the design and management problem would 
include values ranging from 


For these examples, assumed that want apply given depth 
water everywhere the field some desired uniformity. 


Design 


Infiltration time: 160 min. 

Infiltration exponent: 0.5. 

Roughness: 0.15. 

Design uniformity: 80%. 

Available flow rate: 283 L/s (10 cfs). 

Field length: 268 134 (880 440 ft). 

From and 25, calculate and with all units meters and 
seconds: 


LEVEL-BASIN DESIGN 


For field length 268 calculate 268/270 0.993. From Fig. 
12.0 L/s-m (0.129 cfs/ft). The border width 283/12.0 
23.6 (77 ft). The area the basin 0.63 (1.56 acre). 

For field length 134 0.496. From Fig. with 0.496, 
80%, find 0.75. The unit flow rate (0.75)(2.86) 2.15 
L/s/m (0.0234 cfs/ft). The basin width 283/2.15 132 (433 ft). 
The area 1.77 (4.4 acre). 

This example summarized Table halfing the length run for 
this example, the basin area can almost tripled without any sacrifice 
uniformity. Also, decreasing the length run, the distribution uniformity 
can greatly improved for the same basin area, thus potentially decreasing 
the amount water applied. 


TABLE 1.—Design Basin Dimensions 


per 
sec- 
ond- 


Fig. can used estimate the effects application depth for the two 
basin lengths. From and 21, calculate and (from results Table 
with 268 with units seconds and meters 


6/13 
u 


From Fig. this corresponds distribution uniformity 80% (Figs. and 
give identical results). Similar calculations can made for 12.70 and 
7.62 net depths. These results and the results for the shorter basin are 
given Table For the shorter basin, the distribution uniformity does not 
respond drastically changes flow rate. 
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Basin length: 201 (660 ft). 

Basin width: (200 ft). 

Roughness: 0.10. 

Flow rate: 283 L/s (10 cfs) 20%. 


The infiltration water into the soil can change considerably over the season. 
can also change with the soil moisture deficit. For this example, the following 


TABLE 2.—Comparison Basin Designs 


liters 

per 
meter meters 


(1) (2) (6) 


TABLE 3.—Effects Flow Rate and Net Depth Distribution Uniformity and 
Application Time for Level-Basin Management 


liters Net Depth 
per second 


(1) 


relationship assumed for and 10.16 in.), 300 min; 7.62 
in.), 240 min; 5.08 in.), 160 min; all with 0.5. From and 
calculated and distribution uniformities found from Fig. for the three infiltration 
functions and for flow rates 227, 283, and 340 L/s cfs, cfs, cfs). 


min- 
268 12.0 0.15 0.197 73.2 0.645 250 3.42 
0.516 160 2.19 
134 0.15 0.0758 78.7 12.70 1.68 250 3.18 
134 0.15 0.116 47.8 1.09 250 5.23 
(3) (4) 
(a) Distribution Uniformity, percentage 
340 
(b) Application Time, minutes, rounded 
227 102 
283 
340 
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The results are given Table Table shows the uniformities for 227 
L/s, 283 L/s, and 340 L/s for the three infiltration functiors previously given. 
Table gives the resultant application cutoff times. table like this 
can very useful tool for the efficient application irrigation water 
level basins. 

The values distribution uniformity and application time this table are 
estimates for the potential application efficiency. Factors not considered here 
variation the actual values roughness, intake, and flow rate, could 
make these efficiencies unattainable. However, this makes the graphical solution 
presented here even more valuable, since the effects the variability these 
factors can easily evaluated. 

The final subsurface distribution infiltrated water can estimated 
using dimensionless advance curves and volume balance procedure, found 
work the first and second writers (5). Their solution and the one presented 
here differ that the former, recession assumed occur simultaneously 
over the entire field, whereas the curves presented here allows for recession 
occur the tail end the basin first, progressing back the head end. 


The dimensionless solution advance and recession level basins was 
extended show the distribution uniformities for wide variety conditions. 
This was then transformed into two representations distribution uniformity 
that are more useful for designing and managing level basins. One solution 
displays the effects net infiltration depth and the necessary infiltration 


opportunity time distribution uniformity. The other displays the effects 
field length and flow rate distribution uniformity. 
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The following symbols are used this paper: 


basin area; 

exponent for power infiltration function; 
basin width; 

units coefficient for Manning’s 
reference drag per unit plane area; 
distribution uniformity; 

drag per unit plan area; 


A = 
B = 


LEVEL-BASIN DESIGN 


application efficiency; 

miscellaneous, unknown functions; 
dimensionless relative infiltration constant; 
constant for power infiltration function; 
field length; 

Manning roughness coefficient; 

flow rate irrigation stream entering basin; 


flow rate per unit width irrigation stream entering basin; 
field slope; 


reference time; 
time; 
time which irrigation stream cut off; 
dimensionless relative volume water; 
reference horizontal distance; 
horizontal distance along basin; 
advance distance along basin; 
reference surface water depth; 
surface water depth; 
reference depth water infiltrated; 
depth water infiltrated; 
soil moisture deficit; 
minimum net depth water infiltrated; 
infiltration opportunity time for depth z,; and 
miscellaneous, unknown functions. 


~ 


x 


Subscripts 


Superscripts 


indicated dimensionless variable. 
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SIMPLIFIED METHOD FOR RAINFALL INTENSITIES 


INTRODUCTION 


developing countries, rainfall intensity measurements are frequently not 
available, and, even when reported, deficiencies record keeping, reading 
charts and graphs, conversions other units may make the records less 
usable less reliable than desired. good universal relationships can found 
between recorded data and intensity, probability distributions rainfall inten- 
sity-duration-frequency amounts may then found that will useful for 
evaluation records and for many purposes related water management 
planning, design, and research. 

The frequency high-intensity precipitation required for several engineering 
purposes, including agricultural drainage; estimation extreme floods for 
inadequately gaged streams; estimation erosion and reservoir sedimentation; 
and designing the dimensions storm sewers, flood channels, and spillways. 

This study explores the usefulness some rainfall amount-intensity relation- 
ships. Those least partially evaluated include regressions predicting maximum 
monthly rainfall from mean monthly rainfall, 10-yr return period, 24-h duration 
amounts from maximum monthly rainfall, and the full range rainfall intensity- 
duration-frequencies from the 10-yr-24-h values. 

Data from Central America, the United States, India, Australia, and, 
limited degree, the world were used. Additional study and local evaluation 
verification recommended prior use the methods presented, particularly 
areas having rainfall regimes significantly different from those evaluated. 

Maximum Monthly Precipitation, PMX.—The World Meteorological Organiza- 
tion, WMO, (8) gives 30-yr mean rainfall, PM, and 30-yr maximum for monthly 
values, PMX, for many weather stations. Hargreaves (4) used and PMX 
from WMO (8) make world water balance study. Since mean precipitation 
values are published more frequently than 30-yr maximum values, was 
considered worthwhile develop equations for estimating PMX from PM. The 
regression equation the form: 
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TABLE 1.—Regression Coefficients for Estimating 30-yr Maximum Monthly Rainfall 
from Mean Monthly Rainfall 


Regression 
coefficients 


Coefficient 
determination 


Standard Degrees 
deviation freedom DOF 
(5) 


Country 
(1) 
Albania 
Argentina 
Austria 
Bangladesh 
Brazil 
Belgium 
Bulgaria 
Canada 
Chad 
Congo Republic 
Czechoslovakia 
Dahomey 
Denmark 
Dominican 
Republic 
Ecuador 
France 
Germany 
Ghana 
Greece 
Hungary 
Iceland 
India 
Ireland 
Italy 
Ivory Coast 
Japan 
Korea 
Malagasy 
Mali 
Mauritius 
Mexico 
Mozambique 
Netherlands 
New Zealand 
Niger 
Nigeria 
Norway 
Pakistan 
Philippines 
Portugal 
Romania 
Senegal 
Sierra Leone 
South Africa 


* 
241 1.29 145 
1.69 139 104 
112 1.87 176 


RAINFALL INTENSITIES 
TABLE 1.—Continued 


Spain 

Sri Lanka 
Sudan 
Sweden 
Switzerland 
Taiwan 
Tanzania 
Thailand 
Togo 

Turkey 
Uganda 
United Kingdom 
United States 
Upper Volta 
Uruguay 
Yugoslovia 
Zambia 


For each country’s data evaluated, values were recorded and for 
PMX and millimeters well the coefficient determination, 
the standard deviation, and the degrees freedom, DOF. The calculated 
regressions appear Table 

For 69% the countries shown, the regression equation predicts 70% 
more the variance PMX. For those countries that have low values 
R’, assumed that two more distinct types rainfall rainfall distributions 
were probably analyzed together. 

The computer program printed the actual values and the computed values 
PMX for those locations for which there large residual large difference 
between measured and estimated values. The large positive residuals, RPMX, 
were used regression analyses with values PMX estimated from the 
various regression equations shown above. The resultant equation 


RPMX 156 1.14 PMX 


with 95, 130 and DOF 106. The envelope relationship obtained 
adding 180 RPMX. Maximum monthly rainfall will, therefore, very 
rarely exceed 336 1.14 the maximum recorded month during 
30-yr period record. 

Ten-Year, 24-Hour Rainfall, P10,24.—Records daily rainfall are much easier 
obtain than records for depth-durations intensities for shorter time periods. 
Records adequate for estimating 10-yr return periods are more commonly available 
than are longer data series. For purposes this study, the 10-yr return period, 
24-h duration amounts, 10,24 were accepted standard and used for 
developing equation for estimating other return and duration periods. Hersh- 
field (5) shows that 24-h (1,440 min) rainfall averages 1.13 times daily (observation 
day) precipitation. 

regression analysis was made for United States data using values PMX 
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and 10, amounts shown the atlas Hershfield (5). The equation for 
precipitation amounts millimeters is: 


with 84%, and DOF 46. Similar equations with approximately 
equal values were developed computer analysis data from Salvador 
and Nicaragua. Limited data from Honduras indicates that Eq. overestimates 
10,24 for that area. Honduras very mountainous and the relationship appears 
less reliable the lee mountains and areas surrounded mountains. 

daily rainfall records are available, 10,24 can determined probability 
analysis. the recurrence interval, the number years record 
(the number items the annual series) and the rank the items 
the series—m being for the largest daily rainfall, then can estimated 
from one the following equations. (Source Eqs. and WMO No. 
168 80, 1970. Source Eq. ASCE Manual 28, Hydrology Handbook, 


The frequency analyses may carried out with either annual series 
partial series. The annual series uses the largest event each year. The 
partial series uses all values exceeding arbitrary value regardless the year 
occurrence. Hershfield (5) shows linear relationship between annual series 
and partial duration series with partial duration series results equal 1.11 
the annual series. 

Hershfield’s (5) atlas relies the Gumbel procedure for fitting data the 
Fisher-Tippett type distribution determine rainfall amounts for longer return 
periods. This latter procedure depends upon functions the sample size and 
the standard deviation. 

Hershfield (5) indicates that probable errors depth duration rainfall values 
presented the atlas are the order 10-20%. Depth-duration rainfall amounts 
can not easily accurately extrapolated over large areas. Many rainfall 
records are site specific. The average 10-yr event can not determined accurately 
from records. Confidence limits and probable sampling errors are 
estimated assuming particular form frequency distribution. The tabulation 
presented Table showing values used Eq. indicates that the 
skew different rainfall distributions varies greatly. 

would seem logical assume that length record, representativeness 
the site and the skew the distribution typical for that particular climate 
area have more influence the accuracy results than the form the 
assumed distribution. Bell (1) indicates that for 68% confidence limits, one 
standard deviation with record, the upper limit for estimated 50-yr 
period 220 and the lower limit yr. 


1949.) 
n+0.4 5 

7 
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Fourth and Sixth Root Relationships.—In study rainfall intensities Central 
America (Un Metodo Sencillo Para Calcular Intensidades Precipitacion 
Centroamerica Kenneth Vogler and George Hargreaves, paper accepted 
for publication Comite Regional Recursas Hidraulicos, Apartado 718, 
Tegucigalpa, Honduras, 1980, p.) was found that depth duration amounts 
varied approximately the fourth root the time, and the sixth root 
the return period, This partially confirmed equation presented 
Bernard (2) for rainfall intensity, The equation is: 


which coefficient depending for its value locality and 
exponent. Values are presented for that portion the United States east 
longitude 100° The range values 0.73-0.82. Eq. can 


Powell (7) suggests that can approximated from, time hours and 
frequency, The equation is: 


Narayana (6) presents table showing maximum storm rainfall area 
India for durations and days and for various return periods. 
the 10-yr, maximum amount assumed correct, other amounts 
can estimated using the equation: 


which constant, then for return periods 5-100 yr, the values shown 
Narayana (6) can reproduced from Eq. The average difference between 
predicted and published values and the maximum difference 11.2%. 

Eq. was applied the values Hershfield’s (5) atlas for the United 
States. was assumed that the atlas values 10,24 were correct and that 
other amount and frequency values could estimated from the given 10,24 
for the range from 5-100 and from 1-24 Differences between 
the atlas values and estimated values using Eq. averaged 6%. The average 
the regressions was 86%. 

Bell (1) used data from the United States, including Hawaii, Alaska and Puerto 
Rico; Australia; and South Africa, show graphically the relationships between 
10-yr, 1-h rainfall and rainfall for return periods 25, 50, and 
100 yr. each case, the relationship linear and passes through the origin 
linear plot. the rainfall depths for 25, 50, and 100-yr return 
periods are estimated using the estimated values differ average 
2.1% from the graphs presented Bell (1). Data from the United States 
are also used show ratios amounts from various durations one-h rainfall 
amounts. The ratios are follows: 


Ratio 0.29 0.57 0.79 1.00 1.25 
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The 30-min and ratios have average difference from those estimated 
from the values using 5.6%. 

Bell (1) compares depth-duration ratios for the United States, the U.S.S.R. 
and Australia. The ratios are very similar magnitude and range for these 
three large geographic areas. This indicates that ratios may much the same 
for world-wide conditions. 

From the above discussion and evaluation appears that the 10-yr-24-h 
precipitation, 10,24 can reliably determined from annual partial 
series, then errors estimating other duration-return period amounts for durations 
min-4 days and return periods 5-100 will not excessive when 
estimated using Eq. 10. Values estimated are well within the limits accuracy 
attainable, considering limitations sample size, accuracy, and representative- 
ness precipitation data. Eq. has advantage over most other distributions 
because its simplicity and ease computation. 


APPLICATION AND 


World-wide coverage mean monthly rainfall presented Wernstedt 
(8). The data are also available computer magnetic tape, (Climatic Data 
Press, P.O. Box 413, Lemont, Pa. 16851) which makes data manipulation very 
easy. Eq. for example, estimates 30-yr maximum monthly rainfall values 
from mean monthly rainfall. Thirty-year maximum values, for less complete 
coverage, are available from various sources including WMO (3) and Hargreaves 
(4). Eqs. and indicate that the values the 30-yr maximum rainfall for 
locations without suitable records can estimated from mean monthly rainfall. 
method also given for the envelope relationship. Eq. estimates the 10-yr, 
24-h rainfall from the 30-yr maximum, and Eq. estimates the other duration-re- 
turn period values for durations min days and return periods 
5-100 yr. 

developing countries, where rainfall intensity data are inadequate not 
available, the above procedures provide means estimating values that can 
used for planning and design water resource developments. major 
deficiency weakness the method the procedure for estimating 
maximum values from mean monthly rainfall. The coefficients the equation 
and the percentage variation predicted are quite variable. Also, 30-yr record 
may not contain maximum value nearly large average 30-yr return 
period amount may include maximum that equivalent average 
return period several hundred years. comparing records from several 


locations, however, attempt can made either eliminate compensate 
for anomalous storms. 


Summary 


Maximum monthly precipitation can reliably estimated from mean monthly 
rainfall. Outlying and envelope values can also estimated. fairly good 
relationship exists between 10-yr, 24-h rainfall and maximum monthly rainfall 
values. For durations days and return periods 5-100 yr, precipitation 
amounts vary approximately the fourth root the time and the sixth root 
the return period. These relationships apply over wide range climatic 
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regions and possibly are global application. Additional study and statistical 
analyses are recommended order more clearly define reliability, universality, 
and general usefulness the methods presented. 

Care must exercised identifying unusual rainfall events and compensat- 
ing for them the analyses. Short records and inadequate station coverage 
nonrepresentative locations can lead substantial error. Errors from short 
records and nonrepresentative data are frequently sufficient magnitude 
that accuracy estimates will not significantly influenced the type 
distribution used analyzing the data. Additional statistical study recommended 
order determine how well the equations presented herein and other similar 
ones can generally applied world-wide distributions. 
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The following symbols are used this paper: 


regression coefficient; 

regression coefficient; 

coefficient depending locality; 

depth precipitation amount; 

degrees freedom (number records minus 2); 
frequency; 

intensity (amount per unit time); 

constant for each location rainfall station; 

rank item the series for the largest; 
number years record number items the series; 

mean monthly precipitation rainfall; 
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maximum monthly value year series precipitation 
measurements; 

year return period hour duration rainfall; 

percent variance predicted the regression equation; 
standard deviation; 

period concentration (usually min h); and 

return period years. 
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OPTIMAL DESIGN BORDER IRRIGATION SYSTEMS 


INTRODUCTION 


Border irrigation systems, both graded and level, are widely practiced methods 
surface irrigation. Effective designs border systems have frequently been 
based arbitrary constraints and performance criteria. The farmer, the 
owner the farm, interested the highest net benefits from crop production. 
Depending upon the amount water available, the cost production and 
the value the produce, the farmer may may not irrigate the whole farm. 
not sure how much area irrigate obtain maximum benefits. Hence, 
procedure analyze given situation and optimally design the irrigation 
system would facilitate effective on-farm water management. This paper presents 
procedure for optimal design border irrigation systems based maximization 
profit while incorporating system operation constraints and the variables 
the operating system. 


Review 


Hall (4) presented simple procedure optimize the design border irrigation 
systems. But only maximization the application efficiency was considered. 
Vierhout (12) applied differential calculus the optimal design border and 
furrow irrigation systems. Once again, the criteria was maximize application 
efficiency. and Liang (13) presented procedure optimize only the length 
run furrow irrigation system. There are many other variables, such 
inflow rate, time irrigation, and net depth application, which should 
simultaneously included the optimization arrive optimal design 
for the system. 

Marjai (7) developed procedure determine the optimum inflow rate into 
border with the other variables remaining constant. The objective was 
maximize uniformity along the length the border. Karmeli (6) presented 
procedure optimize the irrigation quality parameters, such tailwater ratio, 
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deep percolation ratio, and water requirement efficiency for furrow irrigation. 
Though not difficult extend the same procedure border irrigation, 
the procedure becomes highly tedious and lengthy the number combinations 
the variables increases. Besides, this procedure does not allow for system 
constraints. Similarly, Peri, Hart, and Norum (1979) developed the concept 
optimal irrigation depth, considering the economic loss due deficit and excess 
irrigation along with the distribution uniformity the applied water. This 
procedure does not specify the optimal size the irrigation unit, and the 
operational variables, such the flow rate and the time inflow, achieve 
the desired performance. Recently, Reddy and Clyma (9) presented procedure 
optimize furrow irrigation system design based minimum costs, design 
variables, performance parameters, and system constraints. similar approach 
for borders presented here considering maximization profit after deriving 
the relationship between the design variables and the quality parameters. 


System 


irrigated agriculture maximum profit obtained when the losses are minimum 
and the water requirement efficiency optimum. Water requirement 
efficiency defined percent the amount water made available for 
plant use the water requirement the time irrigation (5). Both can 
obtained simultaneously with proper design and appropriate management the 
system. losses are high, excess costs are incurred providing irrigation 
water. Net benefits from crop production are reduced the water requirement 
met each irrigation less than optimum. Thus, there trade-off between 
maximizing water requirement efficiency and minimizing losses. Therefore, the 
irrigation system should designed for optimum net benefits. 

relationship between yield and design variables must established 
optimize the system design. This was achieved two-step process. First, 
relationship was obtained between water requirement efficiency and the design 
variables using hydraulic model. Second, relationship between yield and 
water requirement efficiency was obtained using crop production model and 
the hydraulic model. Later, these two relationships were combined with 
mathematical programming technique optimize the design freely draining 
graded border and level basin irrigation systems. 

Performance and Design Variables.—Hydraulic simulation the applied water 
important component this optimization model. Conceptually the hydraulic 
model represents accurately the operational conditions the irrigation system. 
Actually the hydraulic model provides the volumes (depth) water that enter 
the root zone, deep percolation, and run off the field. The hydraulic 
model simulates these volumes for given conditions intake, slope, and design 
depth (irrigation requirement) and different combinations the design variables, 
such length run, unit inflow rate, and time irrigation. 

All the system variables can constrained specified limits for given field 
conditions. The constraints were follows: 
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which unit inflow rate into the border, L/s; time inflow 
into the border, min; length irrigation run, maximum 
non-erosive stream size, L/s; Minimum flow rate required, L/s; 
maximum time available per irrigation, minutes; length the field, 
and number lengths run. After defining these limits, the values 
the variables were discretized finite number and simulated using 
the appropriate hydraulic model. The models used this study were those 
Strelkoff and Katopodes (11), and Clemmens and Strelkoff (1). Finer discre- 
tizations increase the cost simulation. The length variable was discretized 
into limited number because there are limited number acceptable alterna- 
tives, i.e., the length run may halved reduced one-third. For each 
combination the variables, the water requirement efficiency, the volume 
runoff, and the deep percolation volume were calculated. Here, only freely 
draining graded borders, and level basins are considered. 

Graded Borders.—Graded borders are well suited soils moderate intake 
characteristics and slopes. Efficient irrigation possible balancing the advance 
and recession water. graded borders, the water frequently freely draining 
the downstream end the field. Thus, runoff water becomes important 
component the irrigation system design. desired, the detrimental effects 
deep percolation may incorporated into the problem. 

After simulating for set given conditions and various combinations 
the design variables—length run, unit inflow rate, and time irrigation—the 
values the quality parameters, such water requirement efficiency, and 
deep percolation and tailwater ratios were obtained. statistical analysis 
the data, the following types relationships were defined with high degree 


correlation, between the quality parameters and the design variables. They 
are 


and 


which water requirement efficiency, percent; tailwater ratio 
the volume runoff divided the total volume applied; deep percolation 
ratio the volume deep percolation divided the total volume applied; 
water requirement the time irrigation, millimeters; and 
are constants that are site dependent. For given conditions and constraints, 
these equations provide the relationships between system performance and design 
variables for graded borders. 

Level Borders.—A border basin defined irrigation unit 
zero slope with the downstream end diked. The tailwater ratio zero because 
there runoff. Here, approach similar the one the previous section 
followed. Thus, the relationship given 
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level border irrigation the water requirement efficiency relationship 
sufficient describe the quality parameters. There runoff. Hence, deep 
percolation the only loss the field. Thus, the deep percolation ratio can 
derived easily from the water requirement efficiency. The derivation 
follows: 


5 u 


which volume deep percolation, cubic meters; and and 
units conversion factors. The deep percolation ratio given 


P 


(K;Q,T7,) 


and provide the relationships between system performance and design 
variables for given conditions and constraints. Now, relationship between 
system performance and yield developed. 

Yield and Performance.—Evaluation the optimum level crop production 
for given irrigation system requires crop production model that defines 
yield function system performance each irrigation during the season. 
The system performance each irrigation each section the field obtained 
from the hydraulic model. The model used simulating the yield presented 
elsewhere (10). The depth water applied was assumed constant for each 
irrigation during the season. Any other sequence depth irrigations may 
specified including empirical experimental approaches defining the depth 
the sequence, both. Once the optimal depth irrigation and the crop 
production model are given, the relative yield the crop function 
different levels water requirement efficiency, constant each irrigation 
during the season, was simulated. Different combinations yield and water 
requirement efficiency were obtained varying the design variables: inflow 
rate, length run, and time irrigation. Other variables such slope, intake 
family, design depth may considered. yield versus water requirement 
efficiency function can developed from the above simulation data. For 
given set field conditions, the relationship relative yield water requirement 
obtained from the simulation shown Fig. The relationship was quadratic 
follows: 


E, 
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which relative yield, and defined 


max 


which actual yield, kg/ha; and potential yield under optimum 


100 


a 


FIG. 1.—Relationship between Water Requirement Efficiency and Relative Yield 
Wheat Crop 


conditions, kg/ha. high coefficient correlation 0.96) was obtained 
between relative yield and the water requirement efficiency. 


Problem formulation important component any optimization. The 
problem defined terms objective function (either minimization 
maximization) and related constraints. For this problem, the objective was 
maximize the profit from crop production, i.e., searching for particular value 
water requirement efficiency and the corresponding values the design 
variables that give the optimum net benefit. The profits were due crop 
production particular field. The costs associated with irrigation system 
design are labor, water, energy, ditch construction, and any negative effects 
runoff and deep percolation. direct costs runoff deep percolation 
can quantified, then their costs are included the increased amounts 
water required. After the cost coefficients and the mathematical relationships 
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the quality parameters are obtained, the problem can formulated shown 
below 


value cost water 
the produce 


cost labor cost cost 
ditch production 
construction 


cost deep 
percolation 


profit coefficient, dollars/hectare; cost water, dollars/liter; 
tare-meter; units conversion factor, 60.0; conversion factor, 1/60.0; 
units conversion factor, .0001; cost deep percolated water, 
fraction the time labor utilized during the 
irrigation time; number irrigations per season; number borders 


the width direction; width the border, width the field, 


deep percolation, function (Q,, D,). 
After substituting the yield-water requirement efficiency relationship, Eq. 


into Eq. 11, and neglecting the cost runoff and deep percolated volumes 
the objective function becomes 


value the produce 
(12) 


cost water cost cost cost 
labor ditch production 
construction 


substituting the relationship between the design variables and the water 
requirement efficiency, Eq. into Eq. 12, the objective function given as: 


value the produce 


cost water cost cost cost 
labor ditch production 
construction 
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and the constraints, are given as: 


n,n,,T, = G, nn, 


which minimum width the border, maximum width 
length the run, and total flow rate available the farm, L/s. 

the above problem, the equality constraints were formulated less than 
equal constraints, conformity with the requirement for the geometric 
programming technique, which explained the next section. The constraint 
formulated that satisfied equality. the solution the problem, 
not assumed that the whole field irrigated. Therefore, the constraint 
the width the field, Eq. 14i, not specified equality. profitable 
the farmer, the solution will satisfy the constraint, Eq. 14i, equality 
and the farmer irrigates the whole field. Though, the less than equal 
constraint could have been specified Eq. 14h, instead Eq. 14i, for the 
present purposes, specified the width constraint, assuming that the 
farmer irrigates the full length the field, thereby adjusting the number 
borders the width direction not profitable irrigate the whole field. 


W, 
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TECHNIQUE 


Generalized geometric programming, which applicable engineering design 
problems this type, most appropriate for the above problem. The same 
technique was presented Reddy and Clyma (9) for furrow irrigation systems. 
Generalized geometric programming, GGP, formulated with objective 
function the form: 


with constraints the form: 


the number variables; and and posynomial functions the 
form: 


i 


i=! 
J 


j=l 


which and terms with positive and negative coefficients, respectively, 
the objective function and the constraints; number terms the objective 
function the constraints with positive coefficients (P,(X)); and number 
terms the objective function the constraints with negative coefficients 
The terms and are defined follows: 


M 
mel 


M 


which and coefficients; and €,,, exponents the variables 
the objective function and the constraints; and x,, system variable. 

Eqs. and are called signomials. signomial defined the difference 
two posynomials. the formulation GGP, the terms with the negative 
coefficients are moved the right the inequality, resulting 
Q,(X). Dividing the inequality the right-hand side term results 
rational function posynomials. The posynomial the denominator reduced 
monomial (single term) process condensation defined Dembo 
(2). posynomial divided monomial still posynomial. This posynomial 
then condensed monomial again condensation. After the monomials 
are obtained, the constraints and the objective function are linearized taking 
the natural logarithm the monomial function. This set equations solved 
linear programming. Convergence the solution the original problem 
obtained additional constraints called the original problem. 
solving the linear program finite number times, optimum solution 
obtained the original nonlinear problem. Being nonlinear programming 
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problem, the global solution cannot guaranteed. Different local optima are 
obtained starting different initial feasible solutions. The maximum all 
the optima considered the global solution the problem. signomial geometric 
programming code has been developed solve the problem. 

The values the variables obtained from the above technique are continuous 
(non-integer). the design irrigation system, some the variables, such 
the number lengths run, sets, and borders the width direction, should 
have integer values. Therefore, different technique was attached the above 
procedure obtain optimal solution terms integers for the above variables. 
The branch-and-bound (3) technique was chosen express the related variables 
integer form. 


System 


The generalized geometric programming technique presented above was applied 
the optimal design border irrigation systems: graded borders and level 


TABLE 1.—Cost Coefficients, System Constants and System Constraints for Freely 
Draining Border 


Parameters 
(1) 
(a) Cost Coefficients 


Cost production, dollars 

Maximum production, cubic hectare 

Cost water, dollars 

Cost ditch construction, dollars /linear meter 


System Constants 


Length the field, 
Width the field, 
Slope the field, m/m 
Roughness the field, 
Depth requirement, 

Number irrigations per season 


(c) System Constraints 


Value 
(2) 
56.8 
494.0 
17.41 
3.00 
40.0 
3.25 
805.0 
402.0 
0.001 
0.024 
76.0 
18.0 
0.2716 
L/s 11.20 
0.92 
158.0 
3,600.0 
67.0 
402.0 
9.0 
m 30.5 
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basins. The two cases are discussed separately example with specific 
conditions. 

Example Border.—The data for this problem presented Table 
Using this data the hydraulic model was simulated for different combinations 
the design variables presented Table relationship the following 


TABLE 2.—Relationship between the Design Variables and the Water Requirement 
Efficiency for Freely Draining Border 


Inflow rate, Time Length Water require- Water requirement 
liters inflow, run, ment depth, efficiency, 
per second minutes meters millimeters percent 


(1) (2) (3) (5) 


form was obtained between water requirement efficiency and the design variables. 
0.07678 0.28299 —0.04829 —0.9938 


very good correlation 0.99) was obtained between the predicted and 
actual water requirement efficiency. comparison actual versus predicted 
water requirement efficiency shown Fig. 


3.28 40.0 91.50 99.50 
3.28 60.00 91.50 76.33 
3.28 80.00 91.50 82.74 
3.28 100.00 91.50 102 65.00 
150.00 91.50 98.22 
4.16 30.00 91.50 96.40 
4.16 60.00 91.50 77.29 
4.16 80.00 91.50 83.58 
4.16 100.00 91.50 102 65.50 
4.62 70.00 183.00 102 59.40 

5.58 60.00 183.00 99.50 
5.58 80.00 183.00 83.55 
5.58 100.00 183.00 88.87 
5.58 150.00 183.00 102 73.10 
7.41 40.00 183.00 69.88 
7.41 60.00 183.00 78.31 
7.41 80.00 183.00 102 62.70 
100.00 183.00 89.90 
7.41 150.00 183.00 99.97 
4.62 80.00 366.00 102 66.70 
4.62 100.00 366.00 83.68 
4.62 120.00 366.00 89.37 
4.62 150.00 366.00 96.05 
4.62 180.00 366.00 99.88 
6.04 70.00 366.00 
6.04 90.00 366.00 144 $6.10 
6.04 120.00 366.00 91.85 

6.04 150.00 366.00 102 72.60 
7.41 100.00 366.00 102 66.00 
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After obtaining the relationships between the design variables and the irrigation 
quality parameters, the problem was formulated terms the cost coefficients, 
system constraints, system constants, and the design variables. the present 
study, the effect deep percolation and tail water are not considered. But, 
when appropriate cost coefficients are available, they can incorporated into 
the optimal design process. For the given situation was calculated 


m 3 


3 


56.8$ 
17.41— 989 dollars /hectare 


which price per unit produce After substituting Eq. 


0.07478__028299, -0.04829 


a 


a 


oO 
— 
o 
WwW 
~ 
o 
a 


Actual efficiency, ,in percent 


FIG. 2.—Actual Versus Predicted Water Requirement Efficiency for Freely Draining 
Border 


and the cost coefficients into Eq. 11, and neglecting the coefficients runoff 
and deep percolation, the objective function becomes: 


value produce cost water 


cost cost ditch cost 
labor construction production 
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substituting Eqs. and (at mm) into Eq. and simplifying, 
the result becomes 


The system constraints are given follows: 


The generalized geometric programming technique was applied the solution 
the above problem. The following optimal values the design variables 
and 13. The water requirement efficiency obtained was 91%. 
The net profit was $9,833 for the field which equals 304 dollars/hectare. 
the optimum, the cost design was 186 dollars/hectare. The application 
efficiency for the optimum was 51%. fact, when the objective maximize 
net benefits, much emphasis cannot given application efficiency. The 
technique presented here selects the optimal water requirement efficiency and 
yield without directly considering the application efficiency. 

defining the constraints, care must exercised selecting proper limits 
for the variables. These limits should not different from the limits used 
simulating the hydraulic model; so, the relationships developed may not 
valid. wide range alternatives for the design variables are considered, 
these same limits should included the hydraulic simulation model. 

Example 2—Level Basin.—The optimization technique was applied for level 
basin. The system cost coefficients, constants, and constraints are given 
Table The hydraulic model was used simulate the relationships between 
the quality parameters and the design variables. The different combinations 
the design variables used the simulation are presented Table For 
the given situation, the following relationships were obtained between the quality 
parameters and the design variables: 


E, = 3,1039°"" T° 9201 
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or E —0.8393 D 
and the equation for given as: 


or R, ak 1 0.8597°'"* D, 50 mm 


P 


good correlation 0.96) between the actual and predicted water requirement 


TABLE 3.—Cost Coefficients, System Constants, and System Constraints for Level 
Basin Irrigation System 


Parameter 
(1) 
(a) Cost Coefficients 
Cost production, dollars 
Maximum production, kilograms /hectare 
Cost ditch construction, dollars /linear meter 
Cost water, meter 
System Constants 
Length the field, 
Width the field, 
Slope the field, m/m 
Roughness the field, 
Depth requirement, 
mm/h* 
mm/h 


(c) System Constraints 


efficiency was obtained. comparison predicted versus actual water require- 
ment efficiency presented Fig. The above relationship was used 
formulate the problem. 

The profit coefficient for the given situation calculated as: 


$2,500/ha 
After substituting the cost coefficients from Table and Eq. into Eq. 


301 
Value 
(2) 
1.00 
1,482.0 
2,500.0 
6.56 
3.0 
40.0 
335 
302 
0.0 
0.15 
28.52 
0.1088 
2.25 
11.20 
0.92 
52.0 
Min 3,600.0 
m 168 
| m 67 
nn m 30.5 
m 9 
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and simplifying, the objective function becomes: 


and the system constraints are given 


w 


TABLE 4.—Relationship between Design Variables and Water Requirement Efficiency 
for Level Basin Irrigation System 


Length Inflow rate, Time Water require- Water require- 
run, liters inflow, ment depth, efficiency, 
meters per second minutes millimeters percent 


(3) 


1.00 50.00 88.50 
1.87 25.00 68.00 
1.87 25.00 52.00 
1.87 25.00 44.50 
1.87 40.00 86.00 
1.00 50.00 114 39.30 
1.87 56.00 102 91.00 
1.87 56.00 114 81.00 
100 1.00 120.00 88.90 

134 3.76 25.00 58.00 

134 3.74 33.00 69.00 

134 3.74 45.00 97.00 

134 4.66 30.00 114 53.00 

134 3.76 50.00 100.00 

168 1.87 120.00 102 74.70 

168 3.74 40.00 102 52.00 

201 2.79 50.00 66.00 

201 3.74 65.00 92.00 

201 4.66 45.00 127 49.00 

201 4.66 50.00 102 66.00 

268 2.79 100.00 80.00 

268 3.74 80.00 87.00 

268 4.66 60.00 79.00 
268 5.58 50.00 68.00 
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applying the generalized geometric programming technique, the following 
optimal values the design variables were obtained. 1.20 L/s; 
efficiency satisfied the optimum was 100%. The maximum net profit the 
optimum was $5,797/per field $573/ha. The application efficiency found 
the optimum was 93%. 

the analysis, the cost runoff and deep percolation water was included 
indirectly the cost water provided excess the requirement the 
root zone. But, the negative effects these parameters (runoff and deep 
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FIG. 3.—Actual Versus Predicted Water Requirement Efficiency for Level Border 


percolation) were not considered because lack appropriate cost coefficients 
for these parameters. The cost coefficient for runoff includes the cost removing 
the excess water from the field, and the negative effects water quality. 
Similarly, the cost coefficient for deep percolation must take into account the 
effects waterlogging and fertilizer leaching crop yield. The cost coefficients 
should given terms dollars /unit volume water. Once these coefficients 
are available, they can incorporated into the optimization process. This would 
help evaluating different management practices controlling the quality 
irrigation return flow while increasing agricultural production. 

implicit assumption the problem was that the farmer can apply the 
total available flow rate one border. consideration was given irrigating 
more than one border time. This was not considered the problem 
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reduce the number variables considered, but can incorporated into the 
optimization process desired. was also assumed that the relationship between 
the yield and water requirement efficiency was identical under graded and level 
basin irrigation systems. 


Summary 


combination approach simulation and mathematical programming was 
used develop optimal system design. First, the irrigation quality parameters 
were developed terms the application system design variables using 
hydraulic simulation model. Second, relationship was developed between crop 
yield and the system quality parameters (water requirement efficiency) using 
crop production function. Maximization net benefits was the objective. 
The value the produce, and the costs labor, water, ditch construction, 
and crop production were considered the objective function. The negative 
effects runoff and deep percolation were not considered. The problem was 
defined terms system variables, cost coefficients and system constraints, 
and the generalized geometric programming technique was applied the optimal 
design border and level basin irrigation systems. The design variables considered 
were the inflow rate, time inflow, length the run, number lengths 
run, width the border and number border widths. The procedure gives 
optimal design under given field conditions. addition, the procedure shows 
the possibility combining simulation and mathematical programming techniques 
optimizing system designs. The technique presented provides guidelines for 
improving existing on-farm irrigation systems for better management the 
scarce resources agricultural production. 
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The following symbols are used this paper: 


infiltration exponent; 
exponents irrigation quality parameter equations; 
cost water, dollars per liter; 
cost labor, dollars per hr; 
cost ditch construction, dollars per linear meter; 
cost production, dollars per hectare; 
cost runoff water, dollars per hectare-meter; 
cost deep percolation, dollars per hectare-meter; 
constant infiltration function; 
water requirement the time irrigation, mm; 
water requirement efficiency percent; 
objective function terms cost coefficients, profit 
coefficient and the system variables; 
system constraints terms system variables; 
number terms objective function the constraints 
with positive coefficients; 
number terms objective function the constraints 
with negative coefficients; 
number constraints problem; 
proportionality constants irrigation quality parameter 
equations; 
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units conversion factor, 0.06; 
units conversion factor, 0.001; 
units conversion factor, 10,000; 
constant infiltration function; 
length irrigation run, 
length field, 

maximum length run, 
minimum length run, 
number variables the problem; 

number irrigations per season; 

number lengths run; 

number borders width direction; 

profit coefficient, dollars 

posynomial functions objective function and the con- 

unit price produce, dollars/kilogram 

total flow rate available farm, L/s; 

unit inflow rate into border, L/s; 

maximum non-erosive stream size into border, L/s; 
minimum flow rate into border, L/s; 

deep percolation ratio; 

tail water ratio; 

correlation coefficient; 

time inflow into the border, minutes; 

maximum time available per irrigation, minutes; 

time variable, minutes; 

terms the objective function and the constraints; 
volume deep percolation, m’; 

width the border, 

width the firld, 
maximum width border, 
minimum width border, 

system variable; 

actual yield, 

relative yield; 

potential yield under optimum conditions, kilograms /hec- 
tare; 

cumulative infiltration, mm; 

fraction the time labor utilized during the irrigation 
time; 

units conversion factor, 60.0; 

units conversion factor, and 

exponents system variables objective function and 
constraints. 
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TESTING PLASTIC FILTER FABRICS 


Yung Hai ASCE, Daryl Simons,’ ASCE 
and Phillip 


INTRODUCTION 


Plastic filter fabrics are made from synthetic fibers. They are two types, 
i.e., either woven fabric featuring distinct open areas, nonwoven fabric 
that consists random fibrous mass uniform thickness with torturous 
paths through the fabric. The fabric fiber consists predominantly polyvinylidene 
chloride, nylon, polyester, polypropylene yarns, filaments, fibers. These 
synthetic yarns give rise the term However, the terms filter fabric 
and cloth filters are synonymous with plastic filter. 

Filter fabrics, being both economical and durable, can used alternative 
graded filter. Filter fabrics are economical because costs both natural 
filter material and construction are increasing. Like granular filters, plastic fabric 
filters are designed highly permeable water and yet constrain soil particles. 

Plastic filter fabrics have been used help protect coastal structures associated 
with rivers, lakes, canals, dams, and drainage systems all types. Principally, 
the filters serve bank protection, subdrainage protection, and means 
protecting the foundation soil. These filters retain their strength undiminished 
long exposure salt fresh water. Many successful applications plastic 
filters have been documented, such the usages plastic filters protect 
beaches Florida (1), control the erosion problems bridges (4) the 
Corps Engineers, and serve road support soft, spongy Alaskan 
soil (7). 

interaction exists between the plastic filter fabric characteristics, soil type 
used, and ground-water conditions. analyzing the effectiveness filter 
fabric, the performance mechanism must investigated considering the afore- 
mentioned variables. establishing conducive testing environment, possible 
study the hydraulic properties plastic filters associated with various soil 
types and water head. the intent this study provide information 
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regarding nonwoven filter fabric performance the hope developing 
adequate set guidelines and criteria for plastic filter fabric selection, particularly 


for subsurface drainage systems where filter fabrics are predominantly the 
nonwoven type. 


Previous Test AND SELECTION 


selecting plastic filter, the fabric permeability should adequate 
pass water from the soil while having pore structure that the same time 
will hold back soil. This rationale suggests that design specifications could 
developed based soil type, ground-water conditions, and filter cloth charac- 
teristics. Although current field applications show plastic filters have promise, 
hydraulic characteristics testing has far failed relate filter and soil 
characteristics with system performance. 

The Corps Engineers one the few agencies document results 
their hydraulic characteristics testing (2). Using soil types known gradation, 
the Corps has established set guidelines pertaining the use plastic 
filters. Their results indicate that filter cloth selection dependent the D,, 
size the soil, the EOS (equivalent opening size), and term denoted 
gradient ratio. Two types tests were conducted the Corps: (1) Filtration 
tests; and (2) gradient tests. 

Based upon analysis the filtration tests, the Corps Engineers recommended 
the following criteria used selecting plastic filter. When the filter fabric 
adjacent granular materials containing 50% less weight silt (material 
little plasticity, passing the No. 200 sieve): (1) The fabric pore size 
estimated EOS should smaller than the the soil where the 
the diameter the particle which 85% weight the soil finer; 
and (2) the open area woven filters should not exceed 50%. 

When the filter fabric adjacent soils having little cohesion containing 
more than 50% silt weight: (1) EOS should larger than the opening 
the United States Standard Sieve No. 70; and (2) open area woven fabrics 
should not exceed 10%. reduce the chance clogging, fabric should 
specified with EOS sieve number larger than the opening United 
States Standard Sieve No. 100 (0.0059 in., 0.15 mm). Filter fabrics require 
sand blanket thickness equal larger than in. placed between 
the filter and the soil for soils with 85% more smaller than the No. 200 
sieve. 

The gradient ratio parameter the Corps Engineers uses filter 
comparison. This parameter developed from permeameter testing and 
result hydraulic gradient manipulation: 


which the gradient ratio; the hydraulic gradient the lowest 
in. soil plus filter cloth; and the hydraulic gradient over the in. 
soil between one and in. above the fabric. Another definition for 
the total hydraulic gradient. the gradient ratio becomes greater than three 
within the Corps Engineers suggests that the filter cloth being plugged. 


(1) 
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The design concept the parameter that plugging occurs, there 
water pressure buildup immediately above the plastic filter. Some filter fabrics 
being tested developed gradient ratio greater than 1.0 when subjected 
silty sand. 

Rosen and Marks (7) studied the behavior plastic filters the University 
Tennessee 1975. Unlike the Corps Engineers’ filtration tests which 
were conducted for their experiment was conducted days 
until the flow through the system became constant. Furthermore, the filtration 
and gradient ratio tests were conducted simultaneously. Twenty soil types were 
tested altogether, however, only one filter fabric was considered, i.e., Mirafi 
140 nonwoven filter). The soil types were produced adding various amounts 
either silt, kaolinite, montmorillonite river sand base soil. 

Testing indicated soil compaction occurred within the soil. This was deduced 
from permeability decreases well analysis hydraulic gradients. The 
percentage soil actually trapped the fabric was quite small and did not 
vary significantly with soil type. The more soil was well graded, the smaller 
the amounts fine particles that passed through the filter cloth. Rosen and 
Marks concluded that the design criteria plastic filters must the same 
those established for conventional aggregate filters. The nonwoven filter 
tested was determined effective subdrainage applications for relatively 
wide range soil conditions. 

Ogink (6) investigated the characteristics the fabric aperture and the blocking 
fabrics sand particles. Willardson and Walker (8) conducted tests evaluate 
the response using filter fabrics drain envelope material problem soils 
needing drainage. These writers, the Corps Engineers, and other investigators 
have provided the public with some basic guidelines filter cloth use. 
still the practice some Corps districts not use filter cloth within critical 
regions. The testing has opened the door for filter cloth applicability; however, 
time alleviate the trial and error methodology applying the filter 
cloths while the same time, broadening the scope the plastic filter fabric 
industry. 


Test Apparatus.—The experimental set-up (5) included large head box, 
small rectangular constant head box, permeameters, manometer board with 
tubes, weighing containers, recirculating pump, and drainage. Fig. shows 
the flow diagram for these apparatus. This system allows for the testing 
filter systems simultaneously under identical hydraulic heads. 

The large head box was constructed intercept the city water and then 
convey this water the small rectangular head box. This large wooden box 
was square with (2.95 sides and (0.61 depth. The purpose 
the large head tank was two fold. The tank was fitted with network 
steamlines that the water could heated and maintained desired warmer 
temperature using thermostat. Also, the tank served stilling basin 
decrease the inflow turbulence and air entrainment. 

The small constant head box, (0.61 width, (1.22 length, 
and (0.61 depth, could moved vertical direction along 
support order vary the hydraulic head. Water supplied each the 
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permeameters via the outlets was maintained constant pressure head 
means the constant head box. Plastic tubing connected the head box 
outlet the permeameters. The cylindrical permeameters were constructed 
clear Lucite 11.5 in. (292 mm) diameter, with three distinct sections 
Fig. The top section consisted gate valve regulate incoming 
water flow and distributor head spread out the flow once the water entered 
the permeameter system. Standpipes were connected the top the permeame- 
ters allow the air escape before any serious accumulations developed. 
The middle section was used contain 4-in. (102-mm) deep soil column. 


Manometers 


Drainage 
Test System 


Filter Cloth & 


Woter & Wire Screen 
Drainage Fine Soil Supporter 
— 


Modified Bottom Section 


b. Permeameter 


FIG. Test System 


The top and middle permeameter sections fit together flush with the plastic 
filter and screen support positioned directly below. third section was attached 
the bottom the screen support close the system from the surroundings. 
order measure the pressure head exerted the soil column and the 
filter fabric, six pressure taps were installed the vertical wall each 
permeameter and connected manometers. The first pressure tap was located 
in. above the soil surface, whereas the second, third, and fourth pressure 
taps lay in. increments within the soil depth starting the soil surface. 
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The fifth pressure tap was located 0.125 in. (3.2 mm) above the filter fabric 
and the sixth pressure tap was secured the bottom section below the filter 
fabric. 

5.0 pump was connected the system recirculate the water. 
was thought that recirculating the water, temperature fluctuations and air 
bubbles would minimized. However, the pump was later removed increase 
the fresh water supplied the system order decrease the amount 
bacterial activity within the system. 

Test Procedure.—In order conduct objective testing the plastic filter 
fabrics, test condition consistency well uniformity must maintained: 
(1) Consistent that all soil types are subject the same initial conditions; 
and (2) uniform that for any one particular test, the soils are equal size 
distribution and density for all permeameters. the two objectives, consistency 
was more difficult maintain. The permeameters are all the same size 
and shape. constant head was applied the system, equally each 
the permeameters. The difficulty lay within the initial soil saturation procedure. 
saturating, water flow was upward direction, thereby pushing the air 
pockets and out the soil sample. Every effort was made keep the 
upward flow slow enough that the soil was saturated uniformly while the 
same time displacing only air pockets. Even so, local piping still occurred. 
testing soil types with fine particulates, was advantageous add 3-in. 
layer pebbles top the 4-in. soil column. This surcharge pebbles 
significantly reduced piping effects. 

The other objective that must considered testing soil uniformity. 
The soil each the permeameters must the same size distribution 
and density. Soil components were mixed small cement mixer for 8-h 
period. This careful mixing soil components which comprised soil type, 
assured approximate uniformity. Once the soil was placed within the per- 
meameters, was compacted appropriate depth that would provide 
nearly maximum density. Three distinct layers the soil column were individually 
compacted. 

The testing procedure began with the heating and chlorination water 
the large head tank. The steam lines within the large head box and the thermostat 
were set that the water temperature within the system was approximately 
room temperature, 68° Once sterilized, lbs (11.3 kg) soil was placed 
each the twelve permeameters, and compacted thickness 4.16 
in. (106 mm). This corresponded soil density 100 (1,600 
The soil was then ready saturated. The permeameters were filled 
upward direction approximately in. (25.4 mm) above the top soil surface. 
The saturation process usually required depending soil composition. 
Once the saturation process was completed, the transport line was connected 
the permeameter inlets. Water entered the system downward direction 
and sprayed against the inside permeameter walls through the distributor heads. 
The inch water above the soil depth during saturation guaranteed dissipation 
the incoming water energy. This eliminated soil surface erosion and other 
soil disturbances. The permeameters were allowed fill and the pressure readings 
stabilize. Once the systems reached stabilized conditions, water inflow was 
adjusted that constant head was applied upon the soil depth each 
permeameter. After adjustments were made, further increment time was 
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necessary for pressure and flow stabilization. When all the permeameters 
indicated the required constant head, testing began. 

The data collection program includes measurements discharges and pres- 
sures. Discharge measurements were collected daily basis from each 
the permeameter outlet taps. Once all discharge and pressure measurements 
were terminated, the permeameter systems were dismantled. The bottom section 
the permeameter served sediment trap for the fine particles which passed 
through the filter fabric. Soil and filter samples were analyzed determine 
the bacterial activity within the soil and filter fabric. The bacteria analysis was 
performed microbiologist and the tested filter fabrics were analyzed with 
scanning electronic microscope determine filter clogging bacterial activity, 
both. 

Tested Soil, Filter Fabrics, and Water Head.—Four soil types were used 
analyze filter fabric performance. These four soil types represent variety 
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FIG. 2.—Size Distributions for Tested Soils 


fine soils which the Corps Engineers felt critical plastic filter 
usage. The soil size distributions are shown Fig. These soils are classified 
sand (Soil Type sandy loam (Soil Type 2), sandy loam (Soil Type 3), 
and loam (Soil Type 4), the Triangular Soil Classification. Kaolinite was 
used the clay component testing. This type clay was chosen for its 
insignificant swelling nature. 

different plastic filter was placed within each filter tube. Table lists 
the various filter fabrics and the corresponding equivalent opening sizes used 
the four soil tests. All the filter fabrics tested were nonwoven fabrics. For 
Soil Type the hydraulic gradient (head difference between pressure Taps 
and divided the soil thickness) imposed the test was 2.5 (762 
mm) water depth upon (305 mm) soil. For Soil Type the hydraulic 
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gradient simulated water depth upon soil. For Soil Types 
and the hydraulic gradient was maintained water depth upon 
soil for 800 After 800 the hydraulic gradient was increased 1/2 


TABLE 1.—Tested Filter Fabrics 


EOS ranges: 

United States Mean opening 
Filter standard size, 
used sieve size millimeters 


(2) (3) (4) 


4 


each day until reached water depth upon soil. Then the 
hydraulic gradient was lowered the original value. 


AND CONSIDERATION 


Method hydraulic analysis the filter fabrics consisted 
six-phase study. These six phases include water filtration discharge study, 


test 
(1) (5) (6) 
80-100 0.17 0.65 1.06 
70-100 0.18 0.69 1.13 
140-170 0.10 0.38 0.63 
40-50 0.36 1.38 2.25 
80-100 0.17 0.65 1.06 
70-100 0.18 0.69 1.13 
40-50 0.36 1.38 2.25 
70-80 0.20 0.77 1.25 
0.42 1.62 2.63 
0.30 1.15 1.88 
100 0.15 0.54 0.88 
80-100 0.17 0.71 17.0 
70-100 0.18 0.75 18.0 
140-170 0.10 0.42 10.0 
40-50 0.36 1.50 36.0 
70-80 0.20 0.83 20.0 
0.42 1.75 42.0 
10A 0.30 1.25 30.0 
3-01 0.95 3.97 95.0 
70-100 0.18 0.72 24.0 
0.42 1.68 56.0 
10A 0.30 1.20 40.0 
3-01 0.95 0.80 126.7 
70-100 0.18 0.78 45.0 
0.42 1.83 105.0 
10A 0.30 1.30 75.0 
3-01 0.95 4.12 347.5 
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water pressure analysis, analysis sediment discharges through the plastic 
filters, soil size distribution analysis, analysis permeability changes, 
and bacteria well fabric fiber analysis. avoid confusion, the following 
terminology used. 


fabric n’’ denotes the tested filter fabric identified the Number 


tube n’’ denotes permeameter device for testing the filter fabric 
identified the number 


Using the measured discharges and pressure distributions, the gradient ratio 
and the permeability coefficients based the average and bottom inch pressure 
gradients were computed. 


The governing equation for soil permeability computation was given Darcy’s 
law 


which the water discharge; the permeability coefficient the 
soil; the hydraulic gradient; and the cross-sectional area the soil. 

The average hydraulic gradient was computed dividing the water head 
applied the filter systems 4.16 in. soil depth. The bottom hydraulic 
gradient was computed dividing the pressure recorded within the bottom 
inch soil in. soil. Knowing the average water velocity and the hydraulic 
gradients the permeability coefficients based upon the average and bottom 
pressure gradients were computed for each filter tube. 

The gradient ratio calculated using Eq. was the index used the Corps 
Engineers indicate that filter fabric was clogging. The value was 
ratio the hydraulic gradient within the bottom inch soil (including the 
filter fabric) the average hydraulic gradient. Due the fact that the average 
hydraulic gradient for the four soils tested was maintained constant value 
during most the testing time, the gradient ratio would therefore indicator 
the pressure buildup within the bottom inch soil. 

Upon dismantlement the filter tubes, core sample the tested soil 
was analyzed determine the bacterial activity within the soil. The tested 
filter fabrics were examined using scanning electronic microscope and analyzed 
microbiologist determine the extent soil plugging and bacterial activity. 
The small amounts fine particulate that passed through the fabrics were 
weighed this time. However, most the soil testing considered large 
concentrations particles that were finer than the opening sizes the plastic 
filters. Therefore much the sediment passing through the fabrics occurred 
when the dry soil was placed the filter tubes. 

Soil Tests.—For Soil tests, was found that each the filter systems 
experiences declining flow rate with time. Fig. plotted using discharge 
data from filter tube shows typical filtration response for Soil tests. 
The discharges declined from the initial value 600 cc/min final value 
approximately cc/min after 1,000 testing. The pressure distribution 
showed definite increase pressure difference after testing between 
the top first and second inch soil. After about brownish-red soil 
layer developed within the top in. soil. When dismantling, this dark layer 
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the soil was found heavy accumulation two bacterial types, some 
mold particles, and small trapped silt and debris. The bacterial types were 
sticky capsular material which could retard the water flow. The laboratory 
environment serves excellent atmosphere for bacterial activity due 
the constant lighting and warm temperatures. However, city water with chlorine 
residual ppm-2 ppm was used testing. Evidently most the residual 
chlorine was dissipated within the constant head tank with the remaining chlorine 
unable microbial activity below the top inch soil. The filter cloths 
had some bacteria attached but under light microscope they did not appear 
damaged deteriorated. scanning electron microscope was used evaluate 
each filter fabric for possible soil bacterial blockage within the fibers. Each 
filter fabric examined showed indication fabric clogging. Samples the 
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FIG. 3.—Typical Flow Rate and Pressure Distribution for Soil and Filter Fabrics 


microscopic photographs are shown Fig. The figure shows the control 
(original fabric) and one filter section both magnified times. The microscopic 
photographs confirmed the findings the light microscope for Soil tests. 

the bacteria grew within the soil, the water velocity was greatly reduced 
indicated the flow rate decrease. This decreased the soil permeability 
accordingly shown Fig. The bottom inch soil experienced smaller 
reduction permeability than the average soil permeability shown Fig. 
believed that the plugging the filter system mainly occurred because 
bacterial activity the top in. soil which caused big pressure drop 
across this soil layer. This shown the reduction gradient ratio shown 
Fig. Initially the gradient ratio was larger than 1.0 due larger compaction 
the lower portion soil. The ratio began decreasing the pressure drop 
increased within the top in. soil. 

Fig. plot the initial and final discharges versus the equivalent opening 
sizes (EOS) the filter fabrics tested. Due data scatter definite correlation 
could initially drawn between the EOS the filter fabrics and the water 
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discharges. equation representing the weighted vertical permeability through 
multiple soils 


FIG. 4.—Photomicrographs Filter Fabric for Soil Test Magnified Times: (a) 
Section; and (b) Test Section 
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soil layer permeability; and and the water head applied the 
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FIG. 5.—Typical Ratio and Permeability Coefficients for Soil and Filter 
Fabrics 
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FIG. 6.—Discharge Versus Equivalent Opening Size for Soil Tests 


soil. The equation shows that the thickness individual layer multiple 
soil system small, the effect this layer upon the measured velocity would 
minimal. this layer was that the plastic filter, initially the filter fabrics 
would have negligible effect upon the filter system. 
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the end testing the discharge response was independent the filter 
opening size. This suggested that the soil (after the growth microorganisms) 
controlled the hydraulic responses the system regardless the type fabric 
used. The sediment amounts that passed through the filter fabrics were less 
than i.e., 0.1 wt% soil permeameter. The mean sediment concentration 
for the Soil tests was below ppm. There was indication that relationship 
existed between EOS and passing sediment weight. 

Soil Tests.—The majority the filter test system was characterized 
nearly constant flow rate about 100 cc/h. Most the pressure differences 
were approximately constant. Fig. shows typical example. This figure shows 
that the pressures measured Tap approached that Tap The two most 
probable causes were that the soil above Tap was disturbed air introduced 
the soil from the manometer lines the soil was disturbed when first saturated. 

The pressure distributions for Soil tests indicate that water head 
were lost within both the top in. and the bottom in. soil. However, 
pressures recorded within the bottom inch were consistently low. believed 
that when this fine soil type was initially placed within the filter tubes some 
fine particles immediately above the plastic filters passed through the fabrics 
leaving layer more permeable soil. When the bottom inch soil was 
more permeable than the other in. soil, the pressure drop within the bottom 
inch would less. 

Based the same reasoning, the permeability coefficient the bottom-inch 
soil was greater than the overall permeability coefficient shown Fig. 
This feature was reflected the changes the gradient ratio with time. The 
initial value gradient ratio was below 1.0. However, the general increasing 
trend gradient ratio coupled with the reduction permeability within the 


bottom inch soil suggests that the lower portion soil was slowly being 
compacted. 
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FIG. 7.—Typical Rate and Pressure Distribution for Soil and Filter Fabrics 
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Microscope examination Soil was found that negligible amount 
bacteria existed within the soil with virtually bacterial attachment the 
filter fabrics. The bacterial activity was controlled Soil tests boosting 
the chlorine concentration the water ppm and initially baking the 
soil 400° kill any viable organisms. Not only was determined that 
the filter fabrics were free bacteria blockage but was also established 
microscopic examination that the fabrics were not clogged soil particles. 

The fact that the water discharge was independent the filter fabric tested 
suggests that Soil controlled the hydraulic responses the filter system. 
was found that Soil Type with large percentage silt and clay was 
less permeable than the filter fabrics tested. Within the tested EOS ranges 
the plastic filters under Soil performed the same. The sediment that passed 


the filter fabrics was not measurable due the small amounts and very fine 
particulate. 
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8.—Typical Gradient Ratio and Permeability Coefficients for Soil and Filter 
Fabrics 


Soils and Tests.—Soils and tests operated the same water head 
Soil tests for the first 800 test period, and then the head increased 
and decreased described earlier. During the constant head test period, the 
flow rates filter tubes decreased from about 250 cc/h for Soil 
tests and from about cc/h-10 cc/h for Soil tests. The reason was attributed 
the bacterial growth the filter fabrics. 

Upon bacterial analysis was found that negligible amount bacteria 
existed within the soil, however significant amount iron bacteria existed 
within the filter fabrics. These bacteria originally propagated upon the metal 
supporter, and then grew into the filter fabric directly above. simulating 
this particular bacterial growth rate special laboratory under ideal conditions, 
was found that the bacteria developed about days. This coincided with 
the reduction the filtration rates. The amount iron bacteria present within 
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the filter fabrics had significant effect the water flow rates. The same 
test preparation procedure applied for Soils and applied for Soil However, 
the metal support screen used beneath the plastic filter was painted with only 
one coat galvanizing compound rather than five coats. This single coating 
galvanizing compound was not sufficient preventing rust accumulation 
upon the metal which then supported microbial growth. 

The water head was increased raising the constant head tank after 800 
testing. direct correlation existed between the water head increase and 
the flow rate. Once the water head was lowered 1150 the original height, 
the flow rate approached constant value approximately equal the discharge 
recorded immediately before the head increase. 

The gradient ratio was initially below 1.0 for the majority the filter tubes 
due disturbance caused backfill saturation before the test. the 
gradient ratio increased above the value 1.0 while the permeability the 
bottom inch the soil decreased from about cm/s-2 cm/s 
for Soil tests and from about cm/s for Soil 
tests. believed that bacteria was developing the filter fabric the 
water velocity reduced and the pressure drop within the bottom inch soil 
increased. Hour 750 the bacteria had significantly affected the discharge. 
After Hour 800 the water head was raised. The permeability coefficients within 
this time range were relatively constant about cm/s for Soil 
tests. This indicates that the velocity through the system directly corresponded 
with the increase hydraulic gradient. Although the filter fabrics accumulated 
bacteria, the soil was virtually unchanged. 

Like the previous two soil tests, correlation between water discharge and 
the fabric EOS could not established for the Soil tests. The aforementioned 
weighted permeability equation shows that the filter fabric had negligible effect 
upon the filter system. Therefore, the soil initially controlled the hydraulic 
responses the system. the end testing, the soil also controlled the 
hydraulic responses. Even though the bacterial growth was initiated the 
metal support screen, the activity spread through the filter fabrics the adjacent 
layer soil. 

Analysis Results.—A filter fabric may serve permeable solid constraint 
real filter. Quite often particular soil layer within multiple layer soil 
system homogeneous soil body, having adequate size distributions, can 
work filter limit the migration fine particles within the soil. this 
case, the soil layer serves soil filter and the filter fabric actually works 
permeable constraint stop adjacent soil particles from washing through 
the fabric. Because there are few fine particles reaching the fabric, the filter 


fabric will not clogged soil. According filtration study, for filter 
fabric constrain large particles 


Filter Pore Size EOS) 
(of soil) 


which P,, pore diameter fabric which 95% finer; EOS the 
equivalent opening size fabric; and the soil particle diameter 
which 85% finer. However, the Corps Engineers suggested that 
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D,, (of soil) 


Our experimental results show that even though the range EOS/D,, varied 
from 0.38 for Soil for Soil test shown Table the loss 
soil particles through filter fabrics was very small. This indicates that Eq. 
too restrictive. Therefore, average ratio proposed: 


EOS 
(of soil) 


average hole size soil about one-fifth the soil particle diameter. 
Therefore, layer soil can stop the passage particle the size 
this particle larger than one-fifth the size particles establishing the 


soil layer. For soil layer become soil filter reasonable condition 
then 


Dso Ds; 


Natural soil usually satisfies this condition. the soil contains more than 30% 


TABLE 2.—Gradations Tested Soils 


Silt and clay 
content, 
(3) (4) 


silt and clay, then the filtration rate and soil migration become small 
that the performance the filter system will not affected filter fabrics. 

all the tests conducted this study was determined that the soil controlled 
the hydraulic responses regardless the type the opening size the filter 
used. This indicates that the four soils tested this study are soil filters. Since 
their size distributions (Table satisfy Eq. this equation considered 
satisfactory criteria for evaluating soil filter. 

When the soil body not filter, then significant amounts fine particles 
can carried water through soil voids. this case, filter fabric works 
real filter limit the loss soil particles through the fabric. avoid 
clogging fabric, necessary allow the passage fine particles. The 
following condition then required: 


Filter Pore Size EOS) 
=> 
(of soil) 


(1) (2) 
1.30 1.08 1.14 
2.00 2.86 4.20 
1.47 1.21 1.87 
3.28 3.15 5.50 
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the discharge large enough carry significant amount fine particles 
and piles them upon the surface filter fabric form cake, this cake could 
reduce water flow locally. 

Comparing the experimental results with the Corps Engineers’ criteria 
described earlier, found that the Corps Engineers’ criteria are too 


conservative. The following procedure recommended for selecting filter 
fabric. 


Determine the soil permeability coefficient from laboratory tests from 
available information. 

Determine the size distribution soil layers protected. 

Apply Eq. determine whether not the soil layer can serve 
filter. 

Select range filter fabric having permeability coefficients equal 
larger than that soil. the soil layer filter, then use Eq. select 
suitable fabric. the soil layer not filter, then use both Eqs. and 
for selecting the fabric. For the latter case, the discharge unusually 
high, then local clogging filter fabric may occur. 


possible that bacterial accumulations may responsible for localized 
plugging disturbed field soil, similar that experienced the Soil test. 
Many subdrainage systems contain perforated iron pipe. has been shown 
Soil Tests and that bacteria which metabolize iron compounds can attach 
the filter fabrics and then penetrate and spread through the soil. Iron bacteria 
has caused drainage problems areas with high iron content. Bacterial problems 
may contribute some failures either fabric graded aggregate filters. 


Plastic filter fabrics should designed let water pass and yet retain soil 
particles. Filter fabrics are relatively inexpensive and are easy install. For 
these reasons they offer possible alternative graded filter system. Because 
long-term field experience still limited, some people are reluctant use these 
new cloths. The United States Army Corps Engineers has used filter fabrics 
more than any other single agency the United States. order evaluate 
the hydraulic performance filters, the Corps Engineers developed their 
own design guidelines 1972. Prior this time, very was known about 
the hydraulic responses the manufactured filter fabrics. Even today, the 
majority engineers who consider using plastic filters refer the design concepts 
developed the Corps Engineers. Work Colorado State University and 
other research centers resulting changes design concepts. 

Four different soil types were used evaluate the filter fabric performances 
within this study. Each soil test consisted measuring the water discharge 
and pressure distribution within number filter tubes. Permeameters containing 
4-in. depth soil plus filter fabric comprised filter system. From the 
measured data, soil permeabilities and gradient ratios were computed with time. 

Based the analysis the laboratory tests and field experience with filter 
fabrics, the following conclusions are drawn. 
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natural soil generally well graded and can serve filter soil 
particles limit the migration fines. Therefore, under this condition, filter 
fabric quite often works actually permeable solid constraint rather than 
filter. this case, clogging fabric fine particles does not occur. 

the permeability the soil less than the permeability the filter 
fabric and the soil filter which restricts the migration fine particles, 
then the soil controls the hydraulic responses the system and different filter 
fabrics perform the same. 

When the soil body not soil filter, then significant amounts fine 
particles can carried water through soil voids. this case, filter fabric 
works real filter limit the loss soil particles through the fabric. 
avoid clogging filter fabric necessary allow the passage fine particles. 
the discharge large enough carry significant amounts fine particles 
and pile them the surface filter fabric form cake, then local clogging 
could occur. 

Bacterial activity within the soil upon the filter can control the hydraulic 
responses filter system. 

For 1,000 continuous testing the filter fabrics tested had detectable 
effect upon the hydraulics the filter system. 

The gradient ratio developed the Corps Engineers indicate filter 
clogging should analyzed for longer period time. The Corps Engineers 
system could experience system instabilities that are especially prevalent 
using finer soils. 

The two basic design criteria prevent soil movement developed the 
Corps Engineers state the D,, the soil larger than the the 
filter, determined the EOS test, when using granular soils containing 50% 
less weight silt and when using soils containing more than 50% 
silt, the P,, should smaller than the opening the United States Standard 
Sieve No. 70. These criteria may apply with coarse sands with large water 
velocities; however, with fine sands containing more than 30% silts nonswelling 
clay, both, the soil will definitely control the system and the standard will 
not apply. Using wide range opening sizes, i.e., EOS greater than 
200, with soils large fine concentrations, very little soil will pass through 
the filter with little fine particle migration within the soil. 

Generally speaking, the Corps Engineers’ criteria too restrictive and 
conservative. design procedure proposed this study which will require 
further verification. 
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The following symbols are used this paper: 


cross-sectional area soil; 

gradient ratio; 

thickness soil layer 

hydraulic gradient lowest in. soil plus filter fabric; 
total hydraulic gradient; 

hydraulic gradient; 

permeability coefficient soil; 

permeability soil layer 

pore diameter fabric which 95% finer; 
water discharge rate; 

measured velocity; and 

waterhead applied the soil. 


d, = 
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DISCUSSION 


Note.—This paper part the Journal the Irrigation and Drainage Division, 
Proceedings the American Society Civil Engineers, Vol. 107, No. IR3, 
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Discussions may submitted any Proceedings paper technical note published any 
Journal any paper presented any Specialty Conference other meeting, the Proceedings 
which have been published ASCE. Discussion paper/technical note open 
anyone who has significant comments questions regarding the content the 
note. Discussions are accepted for period months following the date publication 
paper/technical note and they should sent the Manager Technical and Professional 
Publications, ASCE, 345 East 47th Street, New York, N.Y. 10017. The discussion period may 
extended written request from discusser. 

The original and three copies the Discussion should submitted (220-mm) 
(280-mm) white bond paper, typed double-spaced with wide margins. The length 
Discussion restricted two Journal pages (about four typewritten double-spaced pages 
manuscript including figures and tables); the editors will delete matter extraneous the 
subject under discussion. Discussion over two pages long will returned for shortening. 
All Discussions will reviewed the editors and the Division’s Council’s Publications 
Committees. some cases, Discussions will returned discussers for rewriting, they 
may encouraged submit paper technical note rather than Discussion. 

Standards for Discussions are the same those for Proceedings Papers. Discussion 
subject rejection contains matter readily found elsewhere, advocates special interests, 
carelessly prepared, controverts established fact, purely speculative, introduces personalities, 
foreign the purposes the Society. All Discussions should written the third 
person, and the discusser should use the term writer’’ when referring himself. The 

Discussions have specific format. The title the original paper/technical note appears 
the top the first page with superscript that corresponds footnote indicating the 
month, year, author(s), and number the original paper/technical note. The discusser’s full 
name should indicated below the title (see Discussions herein example) together with 
his ASCE membership grade (if applicable). 

The discusser’s title, company affiliation, and business address should appear the first 
page the manuscript, along with the Proceedings paper number the original 
note, the date and name the Journal which appeared, and the original author’s name. 

Note that the discusser’s identification footnote should follow consecutively from the original 
note. the paper/technical note under discussion contained footnote numbers 
and the first Discussion would begin with footnote and subsequent Discussions would 
continue sequence. 

Figures supplied the discusser should designated letters, starting with This also 
applies separately tables and references. referring figure, table, reference that 
appeared the original paper/technical note use the same number used the original. 

suggested that potential discussers request copy the ASCE Authors’ Guide 
the Publications ASCE for more detailed information preparation and submission 
manuscripts. 
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DISCUSSION 


SEASONAL IRRIGATION APPLICATION 


Errata 


The following corrections should made the original paper: 


Page 227, Eq. 18: The upper and lower limits each integral should 
interchanged. 

Page 229, last sentence: Should read one must choose between 
system with 0.20 (UCC 0.84) and 0.10 (UCC instead 


POLYSTYRENE SPHERES EVAPORATION 


The authors are congratulated their efforts reduce evaporation 
losses from water storage facilities. Their efforts are especially important 
those arid regions where evaporative losses are common, and 
where this loss may greater than the amount stored water beneficially 
used. general, the authors’ presentation was clear and however, 
the writer feels that some clarification needed. 

First, the authors indicate that their polystyrene beads did not spin and, thus, 
did not increase evaporation. They report that this contrast results 
presented Myers and Frasier (10). However, the paper Myers and Frasier 
merely references another report (13), which suggests that floating beads might 
increase evaporation the wind causes the beads spin. Ref. did not 
present any results make any suggestion that their beads were spinning and, 
thus, increasing evaporation. The authors also state that Myers and Frasier 
assumed that the reflection solar radiation the beads was the only factor 
that significantly reduced evaporation. Instead, Myers and Frasier reported that 


1980, William Hart, Donald Norum, and Gideon Peri (Proc. Paper 
15687). 

1980, Tariq Mahmoud, and Natiq Bashi (Proc. Paper 15689). 

Hydro., USDA, Sci. and Education Administration, Agricultural Research, 


Northwest Watershed Research Center, Patti 116, 1175 South Orchard, Boise, 
Idaho 83705. 
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reflection solar radiation was the major factor reducing evaporation. This 
was because the beads became wet, and the exposed wetted area was essentially 
the same that open water surface the same size. Therefore, the 
beads did not act true mechanical barrier, but only slightly restricted water 
movement the surface the cover, and, thus, the reflective properties became 
the major factor evaporation reduction. Beads different sizes and surface 
properties may perform altogether differently. 

Secondly, the evaporation reduction efficiencies presented the authors 
could misleading. One should reminded that the authors’ results were 
obtained using very small evaporation pans under abnormal exposure conditions, 
and for very short time periods. Although 6-month study was mentioned, 
only 3-day 10-day results are shown. Previous studies have shown that 
evaporation reduction efficiencies new material are almost always higher 
than efficiencies obtained for the same materials when exposed for periods 
months (10,14). Therefore, the efficiencies presented the 
authors should used with caution, since long-term values may only 30%-80% 
those presented. 

The size the pans used can also significantly affect results. mentioned 
the authors and demonstrated others (14), wind can have considerable 
influence the amount surface area covered this type material, and, 
thus, the effectiveness the material reducing evaporation. 

order for evaporation reduction studies really meaningful, materials 
should tested tanks ponds similar those being used practice 
store water, and for extended periods time. The length the study should 
determined the expected life and cost the cover material—it could 


require only one storage season several years. 

The writer hopes that the authors will continue their efforts reduce 
evaporation losses valuable water supplies, and will present the results 
their long-term studies future reports. 


Bureau Reclamation, Proceedings, Conference Collaborators Evaporation 
Reduction Research Program, U.S. Dept. Interior, Denver, Colorado, January 30, 
February 1963. 

Cooley, R., and Cluff, B., Pond Evaporation with Perlite 
Journal the Irrigation and Drainage Division, ASCE, Vol. 98, No. IR2, Proc. 
Paper 8953, June, 1972, pp. 255-266. 
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Discussion George ASCE 
and Tilemachos 
ASCE 


The physical problem instantaneous drawdown dealt with the author, 
the basis assumptions described the paper under the heading 
has attracted the interest many investigators the past (5, 
and 11). All the referred investigators have approached the problem 
the basis hydraulic theory, which leads the well-known Boussinesq 
differential equation. This type treatment quite compatible with the 
assumption neglecting the seepage surface. the contrary, the attack 
the problem with the accurate potential theory has take into account the 
seepage surface. 

The assumptions made the author and the restriction small drawdown 
imposed Eq. are the conventional ones for dependable treatment the 
problem with hydraulic theory. The resulting Boussinesq differential equation 
describing the physical problem terms the author’s notation may written 


This nonlinear equation, for small slopes the water table and small variations 
s’, may linearized neglecting the second term the right and replacing 
the flow depth, s’, mean constant value, the first one. Under 


these conditions, the linearized equation may written 


For the instantaneous drawdown case the appropriate initial and boundary 
conditions are 


for 
for 


“September, 1980, David Higgins (Proc. Paper 15696). 

Theoretical and Applied Hydr., Dir. Applied Hydr. Lab., Civ. Engrg. 
Dept., National Technical Univ. Athens, Greece. 

and Research Asst. Theoretical and Applied Hydr., Civ. Engrg. Dept., 
National Technical Univ. Athens, Greece. 
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Using the Boltzman variable 


the solution Eq. with the initial and boundary conditions (43, 44, and 
45) 


s’ x? 1/2 


o 


This solution the same the authors solution for large Eq. 25. 

The detailed study the nature the solution, Eq. 47, interest. For 
with which the water table measurements can made. This means that for 
certain time interval, from the beginning, there certain value the 
aquifer length, within which the drainage process active. This length, 
may defined the active drainage length and determined from 


This dependence active drainage length with time has been tested experimen- 
tally the writers Hele-Shaw model. This model consisted two vertical 
parallel crystal plates long and high, placed about apart, 
with solid horizontal bottom and open the atmosphere the upper end. 
The model was connected, upstream and downstream, two plexiglass tanks. 
Commercial glycerine was used the model liquid. The model and the tanks 
were initially filled with glycerine certain depth, The kinematic viscosity, 
was measured with Engler type viscosimeter. Then very fast drawdown 
the liquid the downstream tank, magnitude s’, was imposed with 
proper arrangements the set-up. The resulting free surface the model 
was observed. More specifically the time, from the beginning the 
drawdown the tank was recorded each station, when infinitesimal 
drawdown started each station. More than runs have been tested. For 
each run, the experimental results have been plotted versus and 
strongly linear dependence between and was found for all runs, 
anticipated Eq. 48. Presumably, the constant proportionality between 
and was different for each run. Finally, the experimental data for several 
runs have been presented the dimensionless form Eq. 49, and are given 

The experimental data confirm Eq. very satisfactorily and, thus, verify 
the interpretation the solution, Eq. 47, given before. For the present experi- 
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ments, was found have value 2.175. 

The first writer and Athanasiou (8) have solved numerically the nonlinear 
Eq. 41, following the suggestions made Polubarinova-Kochina (9) and compared 
these results with the solution the linear one, Eq. 47. For values the 
ratio the differences between the results the two 
solutions are small for engineering purposes. For small and close the 
discharge point, every solution which does not take into consideration the seepage 
surface, always existing the physical problem, will error. 

Verma and Brutsaert (10) attacked the same physical problem, but for finite 
horizontal extent, the aquifer. Using the method finite differences 
they solved the problem the basis accurate potential formulation, taking 
into account the seepage surface, and the basis the nonlinear Boussinesq 
equation neglecting the seepage surface. From the comparison their results 
they concluded that, far the results obtained the nonlinear 


FIG. 6.—Variation Dimensionless Active Drainage Length with Dimensionless Time 


Boussinesq equation give very satisfactory results. This means that for 
4h’ the influence the seepage surface importance. 

The author presented Fig. comparison his experimental results with 
the numerical results obtained his proposed solution expressed Eq. 24. 
The writers extrapolated, also, the whole range presented Fig. 
the curve represented Eq. 25, which the same Eq. 47. interest 
that Eq. fits the experimental data even better than the proposed more accurate 
Eq. 24. reality, Eq. overpredicts drawdown, which quite logical taking 
into consideration that the seepage surface has been neglected the formulation 
the problem, while always present nature and experiments. The 
contradictory result that the proposed more accurate solution, expressed 
Eq. 24, underpredicts drawdown for small 1.0). This underprediction should 
increase the seepage surface could taken into account. Thus, seems 
logical, that while Eq. gives compatible results with experimental evidence, 
the more elaborate treatment, presented the author and ending Eq. 24, 
does not improve the accuracy the analytical treatment. 
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The author states that data from Hele-Shaw models are not appropriate. 
the best knowledge the writers, all analytical treatments similar problems 
have been tested and compared with experimental data from Hele-Shaw models 
(10). 
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